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Notation 

A  —  Plant  matrix 
B  —  Control  matrix 
C  —  Output  matrix 
D  —  Modal  actuator 
G  —  Damping  matrix 

G  -  Optimal  constant  full— state  feedback  matrix 

g  -  Generalized  coordinate 

I  —  Identity  matrix 

J  -  Quadratic  performance  index 

K  -  Suboptimal  gain  matrix,  stiffness  matrix 

M  —  Generalized  mass  matrix 

nA  —  Number  of  actuators 

ng  -  Number  of  sensors 

n  —  Number  of  modes 

m 

Q  -  Modal  weighting  matrix 
R  -  Control  Weighting  matrix 
T  —  Transformation  matrix 
a  -  Control  vector 

V  -  Matrix  of  right  singular  vectors 
U  —  Matrix  of  left  singular  vectors 
x  —  State  vector 

y  —  Output  vector 

V  —  Transformed  control  vector 
T  -  Transformation  matrix 

k  —  Sum  of  K’s 
$  -  Modal  matrix 

^  -  Output  matrix 


VI 


a  -  Singular  value 
S  -  Summation/singular  value  matrix 
£  -  Damping  ratio 

O)  —  Natural  frequencies 


c  -  critical/controlled  modes 
s  —  suppressed  modes 
i  -  index  for  variable 
j  -  index  for  variable 
m  —  index  for  variable 
p  -  position  elements 
v  —  velocity  elements 


*  -  Transformed  matrix 

-+■  -  Pseudo  Inverse 
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Abstract 


'  Direct  output  feedback  control  using  one  or  two  controllers  is  applied  to 
the  NASA  Ground  Test  Facility  offset  antenna  model.  This  is  a  test 
structure  designed  to  have  the  vibration  characteristics  associated  with  large 
space  structures. 

The  control  problem  is  transformed  from  physical  variables  into  modal 
variables  and  reformulated  into  a  first  order  system.  This  system  is 
truncated  to  a  reduced  order  model  with  residual  modes  used  only  in 
performance  evaluation.  Optimal  linear  quadratic  regulator  techniques  are 
used  to  design  the  gain  matrices,  and  full  state  feedback  is  approximated  by 
use  of  generalized  inverses  of  the  observation  matrices.  Spillover  is 
eliminated  through  the  use  of  transformation  matrices. 

The  structure  is  shown  to  be  controllable  with  this  method.  Alternative 
sensor  placement  is  explored,  and  found  to  cause  improvement  in 
performance.  The  torsion  modes  are  found  to  be  particularly  important  to 
the  performance  of  the  structure,  but  need  more  sensing  and  actuation  to  be 
adequately  controlled.  Two  controller  systems  require  more  sensors  and 
actuators  than  available  to  achieve  acceptable  performance  with  this 
structure. 


"  i  i 


CONTROL  OF  A  LARGE  SPACE  STRUCTURE 
USING  DIRECT  OUTPUT  FEEDBACK  AND  MODAL  SUPPRESSION 


I.  Introduction 


Mankind  now  has  the  capability  to  place  large  useful  objects  into  space. 
Seeking  energy,  military  advantage,  financial  gain,  or  to  satisfy  scientific 
curiosity,  activities  in  space  will  continue  to  increase  in  the  future.  Such 
missions  tend  to  require  large  and  flimsy  devices,  with  ambitious  performance 
requirements  that  tend  to  conflict  with  a  need  for  low  vibration  levels. 
When  Large  Space  Structures  (LSS)  such  as  these  are  subjected  to  control 
energy  for  pointing  or  maneuvering,  it  is  difficult  to  avoid  exciting  their 
resonant  frequencies.  The  classical  method  of  dealing  with  this  problem  is 
by  constructing  stiffer  structures  or  by  limiting  the  bandwidth  of  the 
commands  coming  into  the  system.  Since  these  vehicles  will  have  to  be 
lifted  into  space,  it  is  desirable  to  actively  control  vibration  rather  than 
avoiding  the  problem  through  passive  methods  of  material  selection  and 
construction.  From  another  point  of  view,  given  a  LSS,  active  vibration 
control  can  be  used  to  increase  the  performance  level  of  that  vehicle. 

The  LSS  control  problem  is  characterized  by:  Infinite  dimensional  and 
nonlinear  in  theory,  large  dimensional  in  practice,  many  low  resonant 

frequencies,  often  grouped  in  "clumps'1,  passive  damping  is  very  light,  and 
performance  requirements  can  be  very  stringent,  particularly  in  pointing 
applications.  The  LSS  is  represented  by  a  discretized  model  obtained 

through  a  finite  element  analysis.  This  technique  yields  errors  in  the  model, 


particularly  in  the  higher  frequency  modes.  Thus  it  is  undesirable  to  try  to 
control  all  of  the  modelled  modes  even  if  it  were  possible.  One  of  the  major 
concerns  in  LSS  control  is  how  to  best  reduce  the  system  size.  From  a 
synthesis  veiwpoint,  this  means  the  determination  of  which  modes  are  critical 
to  system  performance  and  need  control,  which  modes  need  only  be  protected 
from  the  effects  of  the  controller,  and  which  modes  need  not  be  considered  in 
the  design. 

Many  control  strategies  are  appropriate  for  the  LSS  problem.  Most 
center  on  modal  space,  and  use  constant  gain  optimal  full  state  feedback 
[7,8,9,10].  Since  the  full  state  is  rarely  available  to  the  controller,  observers 
or  estimators  are  used  to  provide  the  missing  information.  The  resulting 
dimensionality  problem  has  been  mitigated  by  use  of  decentralized  control 
[1,9,10,11].  This  technique  uses  multiple  controllers,  each  responsible  for 
controlling  a  group  of  modes  using  lower  order  algorithms  than  would  be 
posible  using  one  controller.  The  controller  structure  is  shown  in  figure  1. 
Each  controller  receives  the  same  information  from  the  sensors.  The 
effectiveness  of  decentralized  cost  optimal  control  has  been  demonstrated  by 
several  researchers  [1,9,10].  They  found  that  transformation  matrices  applied 
to  gain  matrices  eliminates  spillover  and  enhances  controller  performance. 
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Figure  1.  Decentralized  Output  Feedback  Structure 


However,  maintaining  four  states  for  each  mode  raises  concern  for  the 
capability  of  computers  to  handle  the  control  of  such  a  large  structure  on 
line.  This  concern  led  Calico  to  propose  decentralized  direct  output  feedback 
(DOFB).  Thyfault  [10]  developed  this  technique  for  three  or  four  controllers. 
The  advantage  of  DOFB  is  that  it  does  not  require  that  an  estimate  of  the 
state  be  kept  on  line,  so  the  computational  burden  on  each  controller  is 
reduced  or  the  number  of  controlled  modes  can  be  increased.  Since  the 
number  of  modes  that  can  be  controlled  is  always  many  fewer  than  exsist  in 
the  structure,  this  may  be  an  overwhelming  advantage. 


This  thesis  applies  the  method  of  direct  output  feedback  to  the  NASA 
Ground  Test  Facility  test  article.  This  is  a  laboratory  test  article  that  has 
been  designed  to  have  properties  characteristic  of  large  space  structures.  One 
controller  and  two  controller  sytems  were  designed.  This  thesis  seeks  to 
show  that  this  control  method  can  be  successfully  applied  to  this  structure. 
The  object  of  the  control  effort  was  limited  to  increasing  the  damping  ratios 
of  the  controlled  modes  as  much  as  possible  while  maintaining  the  stablility 
of  the  uncontrolled  (residual)  modes.  The  time  response  of  the  system  was 
also  explored. 


n.  LSS  SYSTEM  CONFIGURATION 

The  system  model  used  in  this  study  is  a  finite  element  model  of 
NASA’s  Ground  Facility  for  Large  Space  Structures  Control  Verification  at 
Marshal  Space  Flight  Center.  This  is  a  facility  designed  to  test  LSS  control 
schemes,  and  was  not  designed  to  develop  structures  that  will  actually  be  put 
in  space.  In  order  to  accomplish  these  objectives,  the  test  structure  was 
designed  with  the  pathologies  associated  with  large  space  structure  behavior. 

The  test  structure  consists  of  an  ASTROMAST  beam  suspended  from  a 
set  of  three  gimbals,  the  angular  pointing  system  (APS).  The  APS  is  in 
turn  suspended  from  an  x— y  shaker  table,  the  Base  Excitation  Table  (BET), 
that  can  be  used  to  introduce  disturbances  into  the  system.  In  order  to 
simulate  typical  LSS  behavior,  two  tip  attachments  have  been  built:  an 
antenna  and  a  cruciform.  These  attachments  provide  the  closely  packed 
natural  frequencies  associated  with  typical  LSS  behavior.  Figure  2  shows  the 
structure  with  the  antenna  attached,  as  was  used  in  this  study,  and  figure  3 
shows  the  structure  with  the  cruciform  attached.  Future  plans  call  for  the 
antenna  configuration  to  include  evaluation  of  line  of  sight  control 
capabilities. 

In  the  original  configuration,  the  only  actuators  available  for  control  are 
the  three  torque  motors  associated  with  the  APS.  Under  the  VCOSS  II 
program  (Vibration  Control  of  Space  Structures),  four  additional  actuators 


will  be  added.  These  will  be  proof— mass  x— y  actuators.  Four  additional 
accelerometers  will  be  added  co— located  with  these  new  actuators.  As  it 
stands,  the  apparatus  has  x— y  accelerometers  at  the  BET  and  x— y— z 
accelerometers  at  the  mast  tip.  There  are  three  rate  gyros  attached  to  the 
APS  base  plate  and  three  rate  gyros  at  the  mast  tip.  Thus  after  the 
VCOSS  modifications  the  structure  will  have  15  sensors  and  7  actuators. 
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Figure  3.  Test  Structure  With  Cruciform  Attachment 
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The  sensors  outputs  are  processed  by  a  COSMEC-I  microcomputer, 
which  then  sends  the  data  to  a  Hewlett  Packard  9000  series  microcomputer 
for  control  law  calculations.  The  COSMEC— I  system  can  be  used  to  process 
the  control  law  while  the  HP  collects  the  data  for  poet  processing.  The 
outputs  of  the  accelerometers  are  integrated  in  the  COSMEC-I  system.  The 
control  output  is  sent  by  the  COMSEC- 1  in  analog  form,  from  ±10  volts, 
which  is  the  saturation  level  of  the  torquer  servo  amplifier. 

The  x  and  y  pendulum  gimbals  can  exert  51  N.m  of  torque,  and  the 
azimuth  gimbal  can  exert  10  N.m  of  torque.  The  azimuth  gimbai  was 
included  to  facilitate  various  test  scenarios,  and  because  there  is  significant 
coupling  between  the  torsional  and  bending  modes. 

The  beam  is  a  spare  Voyager  ASTROMAST.  The  astromast  weighs 
approximately  5  pounds  and  is  45  feet  long.  It  is  constructed  of  S— GLASS. 
The  beam  was  launched  in  a  folded  state  and  untwisted  in  space.  This 
design  produces  a  longitudinal  twist  of  about  280  degrees,  which  couples  the 
bending  and  torsional  modes. 

The  tip  package  containing  rate  gyros  and  accelerometers  weighs 
approximately  45  pounds.  This  gives  the  structure  the  characteristics  of 
being  pinned  at  both  ends  rather  than  pinned  at  the  top  and  free  at  the  tip. 

The  analytical  model  was  generated  by  the  Interactive  Structures  Matrix 
Interpretive  System  (ISMIS)  finite  element  program.  The  model  included  28 
nodes  and,  where  appropriate,  each  node  had  6  generalized  coordinates.  The 
position  of  the  modes  on  the  structure  is  shown  in  figure  4.  A  133  mode 
model  was  generated  which  was  truncated  to  a  30  mode  model  for  control 
purposes.  The  modal  frequencies  and  a  characterizaton  of  the  modes  can  be 
found  in  table  1.  Unfortunately,  generalized  mass  and  stiffness  matrices  were 


not  available.  However,  since  the  eigenvectors  were  normalized,  the 
eigenvalues  and  eigenvectors  provide  sufficient  information  to  construct  a 
control  model. 
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Figure  4.  Test  Structure  Node  Placement 


Table  1 


NASA  Test  Model  Modes 

Characterization 


[IIP] 


X  Direction  Free  Body  Mode 
Z  Axis  Free  Body  Rotation 
Y  Direction  Free  Body  Mode 
Torsion 

Y- Plane  Pendulum 

X- Plane  Pendulum 

Y- Plane  bending 

Bending  X— Plane  +  Torsion 

Y- Plane  bending 

Bending  X— Plane  4-  Torsion 

Dish  Rim  Motion 

Dish  Rim  Motion 

Bending  X- Plane  +  Torsion 

Y- Plane  bending 

Bending  X- Plane  +  Torsion 

Dish  Rim  Motion 

Dish  Rim  Motion 

Y- Plane  bending 

Local  CW  Arms 

Local  CW  Arms 

Bending  X- Plane  +  Torsion 

Y- Plane  bending 

Bending  X- Plane  +  Torsion 

Local  Antenna 

Bending  X- Plane  +  Torsion 
Y- Plane  Bending 
Y- Plane  Bending 
Bending  X- Plane  +  Torsion 
Dish  Rim  Motion 
Dish  Rim  Motion 


I 


The  main  controller  performance  criteria  was  the  increase  of  damping  as  well 
as  the  maintenance  of  stability.  Time  response  was  also  studied,  but  not 
used  in  the  controller  design  process.  The  time  response  of  the  antenna 
attachment  point,  node  14,  was  used  as  representative  of  structure  time 
response.  The  control  method  used  is  presented  in  the  next  chapter. 


III.  SYSTEM  MODEL 


EQUATIONS  of  MOTION 

Generally,  the  equations  of  motion  for  a  large  space  structure  are  given 


as 

Mg  +  Eg  +  Kg  =  Da  (1) 

where  M  is  a  nxn  symmetric  mass  matrix,  E  is  a  nxn  symmetric  damping 
matrix,  K  is  a  nxn  symmetric  stiffness  matrix,  D  is  a  nxm  matrix  of 
actuator  locations,  f  is  a  nxl  generalized  coordinate  vector,  and  u  is  a  mxl 
control  input  vector.  Equation  1  represents  a  second  order 
eigenvalue— eigenvector  problem.  The  nxn  matrix  of  eigenvectors,  $,  is 
normalized  with  respect  to  the  mass  matrix,  M,  such  that 

4>TM<t>  =  I  (2) 


which  can  be  shown  to  yield  [6] 


4>TE4>  = 


\ 

\ 


$TK$ 


(3) 


where 


I  =  nxn  identity  matrix 
2$<j)  =  nxn  diagonal  damping  matrix 
Ci)2  =  nxn  diagonal  matrix  of  eigenvalues  of  Eq  (1) 


Defining  the  modal  coordinates,  f>,  by 


§  = 


(4) 


Eq  (1)  can  be  transformed  to  modal  coordinates  by  pre-multiplying  by  <f>T, 
yielding 

rj  +  Z^Oirj  +  C02fj  =  $TDq  (5) 

Define  <t>i  as 

<f>a  =  $TD  (6) 


This  can  be  put  into  state  space  representation  of  the  form 

2  =  Ax  +  By  (7) 

where  A  is  the  2nx2n  plant  matrix,  B  is  the  2nxm  input  matrix,  x  is  the 
2nxl  state  vector,  and  a  is  the  mxl  control  vector.  These  are  related  to  Eq 


The  output  equation  is  given  by 

y  =  *j  f  'frpg  =  F  (9) 

In  this  study,  only  rate  sensors  were  considered,  so  the  output  equation 
is  given  by 

y  =  =  4»sf?  =  [  0  |  (10) 

In  modal  coordinates,  the  A  matrix  contains  no  coupling  between  modes. 
The  only  coupling  is  in  the  actuator  and  sensor  matrices.  Thus  the  set  of 
system  model  equations  can  be  split  into  multiple  sets  of  equations,  one  for 
the  controller,  one  for  the  suppressed  modes  or,  in  a  two  controller  system, 
the  second  controller,  and  one  for  the  residual  modes.  A  controller  has  the 
dual  responsibility  of  controlling  the  modes  assigned  to  it  and  "protecting"  its 
modes  from  the  effects  of  observation  or  control  spillover. 

Modal  Control 

The  advantage  of  modal  control  is  that  before  control  there  is  no 
interaction  between  modes.  The  only  linkage  is  through  the  output  equation 
and  through  the  applied  control.  This  means  that  an  n  mode  system  could 
be  arbitrarily  partitioned  into  from  1  to  n  uncoupled  second  order  differential 
equations.  For  the  purposes  of  this  study,  the  system  equation  is  partitioned 
into  three  sets  of  modes.  These  are  the  controller  l  modes,  the  suppressed 
modes  or  the  controller  2  modes,  and  the  residual  modes. 

Thus  the  system  model  can  be  written  as 


£c  =  Acxc  +  Bcu 

x,  =  Afx,  +  Bsu 
fcr  =  A  x.  +  B  u 


(11) 


for  which  the  output  equation  is 


y  =  Ccxc  +  Csxs  +  Crxr 


(12) 


Thus  it  can  be  seen  that  the  output  equation  cannot  be  easily  decoupled  and 
that  the  three  differential  equations  that  make  up  the  suystem  are  coupled 
by  a  common  control  vector,  q. 

The  control  method  employed  in  this  study  is  direct  output  feedback, 
which  feeds  back  the  sensor  outputs  multiplied  by  a  constant  gain  matrix  as 
the  control  input.  Thus  the  controller  gain  must  deal  with  the  coupling  of 
the  modes  by  the  control  (control  spillover)  or  by  the  sensors  (observation 
spillover).  It  will  be  shown  later  that  one  controller  need  not  eliminate  both 
forms  of  spillover.  If  multiple  controllers  are  used,  one  controller  will 
eliminate  control  spillover,  one  controller  will  eliminate  observation  spillover, 
and  any  other  controllers  will  eliminate  both  types  of  spillover.  The  next 
part  of  this  chapter  will  deal  with  the  gain  selection  method.  The  approach 
taken  in  this  study  is  suboptimal  control  using  direct  output  feedback.  The 
method  for  selecting  the  gains  will  be  discussed  in  the  next  section. 

Suboptimal  Control 

Optimal  control  theory  is  used  to  select  the  controller  feedback  gams. 
Linear  Quadratic  cost  optimal  control  (LQ  control)  has  many  advantages  in 
the  synthesis  of  multiple- input,  multiple  output  control  laws  which  suggest 
LQ  techniques  for  the  LSS  control  problem.  However,  LQ  control  requires 
full  state  feedback.  Since  direct  output  feedback  doesn’t  provide  access  to  all 
of  the  states,  we  cannot  use  the  optimal  gain  directly.  Thus  it  is  necessary 
to  operate  on  the  optimal  gain  matrix  to  get  an  appropriate  suboptimal  gain. 


To  obtain  the  feedback  gain  matrix  G,  a  quadratic  cost  function  is 


defined  as 

no 

J  —  1/2 J’ 0  xJQxc  +  aTRodt  (13) 

where 

Q  is  an  n  x  n  positive  semidefinete  weighting  matrix 
R  is  an  m  x  m  positive  definite  weighting  matrix 
and  xc  is  constrained  by  Eq  11a.  For  this  cost  functional  the  optimum 


feedback  gain  is 

G  =  -R-'BJS  (14) 

where  S  is  the  solution  to  the  algebraic  Riccati  equation 

SAr  +  AjS  -  SBCR-'BJS  +  Q  =  0  (15) 

Thusly  the  applied  optimal  control  would  be 

0  -  Gxc  (16) 

yielding  the  closed  loop  system 

*c  =  (Ac  +  BCG)XC  (17) 


This  solution  requires  full  state  feedback  for  proper  application.  Since 
the  output  y  cannot  provide  the  full  state,  an  approximation  to  this  optimal 
gain  must  be  made  for  use  with  direct  output  feedback.  Equating  the  closed 
loop  system  equations  for  direct  output  feedback  and  full  state  feedback 


(Ac  +  BcKCc)xc  =  ( Ac  +  BcG)xc 

(18) 

Ac  +  BcKCf  =  Ac  +  B;G 

(19) 

B,KCr  =  BCG 

(20) 

G 

n 

il 

O 

(21) 
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If  the  Cc  output  matrix  is  full  rank  (making  C.CJ  invertable),  the  suboptimal 
gains  can  be  solved  for  directly 


KCcCj  =  GCj 

(22) 

K(CfCj)(CcCcT)-‘  =  GCj(CcCcT)-> 

(23) 

K  =  GCj(CcCcT)-‘ 

(24) 

When  Cc  is  not  full  rank,  K  is  found  using  the  inverse  of  the  singular  value 
decomposition  of  Cc,  C+,  which  is  called  the  generalized  or  Penrose  inverse. 
From  (21) 

KCcC+  =  GCc+  (25) 

K  =  GCc+  (26) 

Using  one  of  these  methods,  the  suboptimal  gains  can  be  determined  for  an 
output  matrix  of  any  rank. 

Direct  Output  Feedback  Control 

The  system  as  now  equipped  can  most  reasonably  be  controlled  using  a 
one  or  two  controller  model.  Then  there  are  four  classes  of  modes: 
coutrolled/critical  inodes,  suppressed  or  second  controller  modes,  residual 
inodes,  unmodelled  residual  modes.  The  residual  modes  are  used  in  the 
truth  model  if  at  all.  Suppressed  modes  are  not  actively  controlled,  but  the 
feedback  gains  are  selected  to  reduce  observation  and  control  spillover 
between  these  modes  and  the  controlled  modes.  In  a  two  controller  system, 
the  feedback  gains  are  selected  to  cause  minimal  interaction  between  modes. 


The  system  equation  of  motion  is  rearranged  as  follows 

A,  0  0  x,  B, 

0  A2  0  x2  f  B2  |a)  (27) 

0  0  Ar  xr  Br 

The  output  equation  is 

*  =  [  Ci  C2  Cr  ]  [ 

*2  (28) 

*r 

The  proposed  control  for  one  controller  is  of  the  form 

u  =  Ky  (29a) 

For  two  controllers,  this  becomes 

u  =  icy  (29b) 

Where  /c  =  K,  +  K2.  K,  and  K2  axe  the  direct  output  feedback  gain  matrires 
for  system  l  and  2,  respectively.  Thus  the  closed  loop  state  equation  is 

A,  fB,KC,  BtKC2  B,KCr 
B2KC!  A2 1 B2KC2  B2KCr 
BrKC,  BrKC2  ArFBrKCr 

It  is  obvious  that  the  above  control  method  would  not  guarantee  stability, 
even  without  the  residuals.  Since  direct  output  feedback  is  suboptimal,  there 
are  no  guarantees  of  stability,  but  stability  can  normally  be  attained  through 
the  gain  selection  process.  Although  it  is  not  possible  to  avoid  changing  the 
eigenvalues  of  the  residual  modes,  stability  can  be  maintained  by  proper 


selection  of  modal  weightings  and  proper  selection  of  which  modes  to  treat  as 
residual.  To  guarantee  that  the  eigenvalues  of  system  1  and  2  are  as 
designed  (ignoring  the  residuals),  the  gain  matrix  K  or  K  must  meet  a 
number  of  conditions.  The  method  by  which  Eq.  30  is  put  in  the  proper 
form  is  detailed  in  the  next  section. 

Transformed  DOFB  Control 

As  was  stated  in  the  last  section,  the  feedback  gains  must  be  chosen  such 
that  Eq  30  is  upper  triangular.  This  can  be  acheived  by  eliminating  control 
and  observation  spillover  between  the  controlled  and  suppressed  modes  (for 
multiple  controllers,  the  suppressed  modes  of  one  controller  are  the  controlled 
modes  in  the  other  controllers).  For  one  controller  this  means  that  the  1-2 
and  2-1  partitions  of  Eq  (direct  4)  be  zero,  or 
KC2  =  0 

and  (31) 

KCe  #0 

For  two  controllers 

KtC2  =  0 
=  0 

while  (32) 

^  0 

KjC,  0 

While  it  would  be  desirable  to  force  KjC,.  =  0  and  BrKs  =  0,  this  is  an 
impossibility  since  these  modes  are  assumed  to  be  unmodelled.  They  are 


present  here  to  provide  a  system  truth  model. 

From  Eq  (30),  it  is  obvious  that  any  control  that  ignores  the  residual 
modes  cannot  guarantee  overall  system  stability,  since  the  spillover  may 
cause  instability.  This  is  true  of  all  methods,  since  by  nature  of  the  problem 
(infinite  dimensionality)  there  will  always  be  unmodelled  modes.  The  only 
approach  to  this  is  to  try  to  avoid  expending  any  high  frequency  control 
energy  and  to  include  the  modelled  residual  modes  in  any  performance 
evaluation. 

The  method  of  acheiving  the  conditions  given  in  Eq  (31)  will  be  to  find 

transformation  matrices  T  and  T  and  define  a  transformed  gain  matrix 

k,  =  K;r 

and  (33) 

K2  =  tk; 

Then  the  feedback  matrix  for  a  single  controller  is  given  by  K(  and  the 
feedback  matrix  for  a  two  controller  matrix  is  given  by  K  =  K,  + 

The  K*  matrices  are  given  by 

K,*  =  G,(rc,)+ 

and  (34) 

k;  =  g-c  t 

G2  is  the  Riccatti  solution  when  B2  =  TB2  is  used  instead  of  B2.  The 

modifications  to  the  gain  matrices  are  necessary  to  make  the  gain  and 

transformation  matrices  conform.  This  is  discussed  further  in  the  next 


For  one  controller  the  dosed  loop  system  equation  is  given  by 


■  ' 

A,  t  BrKCr 

0 

BrKCr 

*r. 

K 

-- 

B,KCc 

A, 

B.KCr 

L 

BrKCc 

0 

Ar+BrKCr 

■ 

two  controller  system  the  state  equation  is 

given  by 

0 

B^'C, 

- 

*c 

*2 

= 

B^C,  A2fB2/c’C2  B2AC*Cr 

x2 

K 

Br/c*Ct  Br/c*C2 

Ar+BrKCr 

*r . 

(35) 


(36) 


The  sensor  and  actuator  requirements  to  acheive  the  above  closed  loop 
systems  are  given  in  the  next  section  of  this  chapter.  The  transformation 
technique  is  outlined  in  the  next  chapter. 


Sensor/ Actuator  Requirements 


To  satisfy  the  conditions  of  Eq  (31),  the  column  space  of  KT  must  be  in 
the  null  space  of  Cj.  The  effect  of  suppression  is  that  the  eigenvalues  of  the 
suppressed  modes  are  unchanged.  In  a  one  controller  system  this  can  be 
acheived  by  suppressing  either  control  spillover,  where  BSK  =  0,  or  by 
suppressing  observation  spillover,  where  KCS  =  0.  It  is  possible  to  have  two 
sets  of  suppressed  modes.  Then  K  would  have  the  property  that  it 
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eliminated  control  spillover  into  one  set  and  it  eliminated  observation 
spillover  into  the  other  set.  It  is  also  possible  to  eliminate  control  and 
observation  spillover  into  one  set  of  modes.  This  means  that  these  modes 
would  not  be  excited  by  the  controller.  However,  for  this  study  it  made  the 
most  sense  to  suppress  observation  spillover  into  the  suppressed  modes, 
because  there  are  many  more  sensors  available  than  actuators.  In  a  three  or 
more  controller  system,  it  is  necessary  to  eliminate  both  control  and 
observation  spillover  for  some  of  the  controllers. 

The  number  of  modes  that  the  controller  suppresses  is  directly  related  to 
the  number  of  sensors  available.  There  must  be  more  sensors  than 
suppressed  modes  using  the  suppression  method  outlined  below.  Another 
concern  is  that  the  observability  of  the  controlled  modes  goes  down  as  more 
modes  are  suppressed.  This  is  not  demonstrated  rigorously  here,  but  trial 
and  error  computer  runs  demonstrate  the  point.  Direct  output  feedback  has 
more  problems  than  observer  methods  in  this  regard.  Similarly,  using 
actuators  to  suppress  modes  reduces  the  controlability  of  the  controlled 
modes. 

If  there  are  n,  suppressed  modes  and  nsen  sensors,  then  Cg  will  be  n5eoxns. 
However,  for  velocity  feedback,  Cs  =  [  0  |  CVf  J,  so  the  only  condition  needed 
for  suppression  is  that  KCVI  =  0.  The  rank  of  K  cannot  be  greater  than 
the  dimension  of  the  null  space  of  Cvg.  The  null  space  of  CVJ  has  dimension 
n^H  —  r  where  r  <  n,.  Thus  the  designer  of  such  a  controller  would  choose 
to  suppress  fewer  modes  than  he  has  sensors. 

As  has  been  mentioned,  the  number  of  modes  that  are  suppressed  has  an 
effect  on  the  observability  of  the  system.  Thus  it  is  usually  good  practice  to 
suppress  only  those  modes  that  are  affected  significantly  by  the  controller. 


IV.  Suppression  Technique 

The  method  used  to  transform  the  gain  matrix  such  that  it  is  orthogonal 
to  Ca  is  the  same  one  as  used  by  Aldrige,  Miller,  Thyfault  and  others 
[1,9,10,11].  The  optimal  gains  are  modified  such  that  the  closed  loop  system, 
Eq  (30),  yields  Eq  (31),  which  is  in  block  triangular  form.  The  modification 
is  carried  out  by  a  linear  transformation  that  maps  the  gain  matrix  K  into 
the  null  space  of  Ca.  This  linear  transformation  has  been  called  T  in  the 
previous  development. 

T  rids  the  system  of  observation  spillover  between  the  suppressed  and 
controlled  modes.  It  is  also  possible  to  force  the  columns  space  of  K  to  be 
in  the  null  space  of  Ba,  which  eliminates  control  spillover  between  the 
suppressed  and  controlled  modes.  The  transformation  matrix  that  suppresses 
control  spillover  has  been  called  T  in  the  previous  discussion. 

In  a  one  controller  system,  either  type  of  suppression  has  the  effect  that, 
for  no  residuals,  the  eigenvalues  of  the  suppressed  modes  are  unchanged  from 
their  uncontrolled  values.  Suppression  with  one  controller  requires  only  one 
transformation  matrix.  In  contrast,  multiple  controller  models  require 
suppression  of  both  control  and  observation  spillover. 

For  the  suppression  of  a  set  of  modes  while  controlling  another  set  of  modes, 
the  transformation  matrix,  T  or  T,  must  meet  the  conditions 

rc,  =  o 

and  (37) 

rct  *  0 


2 


For  two  controllers  this  becomes 

rc2  =  o 

and  (38) 

TC,  *  0 
while 
BjT  =  0 

and  (39) 

B2T  *0 

As  was  previously  discussed,  there  will  always  be  residual  modes  that 
cannot  be  suppressed,  since  it  is  impossible  to  have  an  infinite  number  of 
sensors  or  actuators.  It  is  also  likely  that  a  large  structure  will  have 
noncritical  modes  that  need  not  be  controlled  to  meet  the  control  objectives. 
In  any  case,  there  are  limits  on  the  number  of  actuators  and  sensors  in  any 
problem  of  interest.  The  effect  of  the  controller  on  the  residual  modes  can 
be  minimized  by  limiting  the  bandwidth  of  the  controller  and  by  reducing 
the  control  energy  that  is  fed  into  the  higher  frequency  controlled  modes. 
Other  studies  [1,2,9,10]  have  considered  multiple  controller  models  where  all 
modes  were  either  controlled  or  residual.  Thus  control  energy  was  applied  to 
modes  very  close  in  frequency  to  the  residual  modes.  One  controller  models 
have  the  property  that  the  suppressed  modes  form  a  barrier  between  the 
controlled  modes  and  the  residual  modes  where  very  little  control  energy  is 
expended. 

The  gist  of  this  discussion  is  that  a  means  of  finding  a  basis  for  the  null 
space  of  B,  or  Cf  is  necessary  to  find  the  T  and  V  matrices.  One 
convenient  source  of  these  vectors  is  provided  by  the  technique  of  singular 
value  decomposition.  Briefly  stated  [3],  for  all  mxn  matrices  C,  there  exist 


unitary  matrices  U  and  V  such  that 
2  =  UTCV 


where  (or  tn>n 


Where  of  are  the  eigenvalues  of  CTC.  The  positive  square  roots  of  of  are 
called  the  singular  values  of  C.  U  is  the  mxm  matrix  of  eigenvectors  of 
CCT  and  V  is  the  nxn  matrix  of  eigenvectors  of  CTC.  If  C  is  of  rank  k, 
then  there  will  be  k  non-zero  singular  values  of  C.  The  bottom  m-k  rows 
of  U  map  the  columns  of  C  to  zero.  These  m— k  rows  of  U  can  be  used  to 
suppress  the  C  matrix. 

Normally,  C  will  be  of  full  rank  so  k=nIup.  C  will  be  less  than  full  rank 
if  its  columns  are  not  linearly  independent.  In  this  case  k  <  n(up.  This 
corresponds  to  redundant  sensors. 

The  above  development  shows  that  T  =  [uf+1  ...  u£]T  can  be  used  to 
suppress  the  C  matrix  given  that  the  system  has  more  sensors  than 
suppressed  modes.  However,  since  T  is  (ns,n-n,up)xnMn,  this  means  that  the 
K  matrix  must  be  nilctx(nsen-nJup)  in  order  to  conform  with  I\  The  gain 
matrix  can  be  modified  to  meet  this  requirement  as  will  be  shown  later  in 


this  section. 


As  has  been  stated  before,  optimal  control  using  a  linear  quadratic 
regulator  requires  full  state  feedback.  However,  direct  output  feedback  gets 
around  this  requirement  by  using  a  pseudo  inverse  as  was  discussed  in  the 
last  chapter.  For  a  one  controller  model,  the  calculation  of  the  gain  matrix 
K*  is  oaly  slightly  different  than  that  of  the  untransformed  gain  matrix  K. 
If  C*  is  defined  as 

c;  =  rcc  (4D 

it  can  be  seen  that 

K*  =  GC*+  (42) 

where  the  optimal  gain  matrix  G  is  the  same  as  given  in  Eq  (14),  which 
gives 

G  =  -R-'BCTS  (43) 

S  is  the  solution  to  the  algebraic  Riccati  equation 

SA  +  ATS  -  SBR“lBTS  +  Q  =  0  (44) 

Since  Eq  (44)  does  not  involve  the  C  matrix,  which  is  all  the  transformation 
technique  changes  for  one  controller,  it  is  reasonable  to  use  the  same  G 
matrix  as  before  transformation.  With  multiple  controllers  it  is  necessary  to 
recalculate  the  G  matrix  [10]. 

Then,  as  in  the  previous  chapter 


(A  +  BK*C')x  =  (A  +  BG)x 

(45) 

A  +  BK*C*  =  A  +  BG 

(46) 

BK*C*  =  BG 

(47) 

K*C*  =  G 

(48) 

Which  yields  Eq  (42). 
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V.  Computer  Analysis. 

A  modified  form  of  the  ACOSS  program  [1,9,10,11]  was  used  to  design  the 
controller.  A  data  flow  chart  of  this  program  can  be  found  on  the  next 


page,  and  a  listing  can  be  found  in  Appendix  A.  After  the  system  model  is 
set  up  and  partitioned,  the  optimal  full  state  feedback  gain  matrix  is 
calculated.  As  an  aid  to  controller  design,  the  loop  is  closed  with  this  gain 
(with  no  suppression). 

The  suppressed  modes  C  matrix  is  then  input  to  an  International 
Mathmatical  and  Statistical  Library  (IMSL)  routine  that  finds  its  singular 
value  decomposition.  This  routine  returns  the  gamma  matrix.  The  gamma 
matrix  is  then  multiplied  by  the  controlled  modes  C  matrix,  and  the 
resulting  matrix,  C*,  is  sent  to  the  gain  Hading  routine.  The  new  gain  is 
multiplied  by  T  to  form  the  transformed  gain  K*.  The  control  loop  is  closed 
with  this  gain  and  the  results  are  sent  to  an  eigenvalue  analysis  routine. 
Finally,  the  time  response  is  optionally  calculated. 


Program  Data  Flow  Diagram 


VI.  Investigation 

The  investigation  mainly  proceded  in  pursuit  of  a  stable  system.  The 
first  cut  at  controlling  the  structure  consisted  of  putting  the  first  11  bending 
modes  in  the  controller,  suppressing  the  next  8  modes,  and  treating  the 

remaining  8  bending  modes  and  the  free  body  modes  as  residuals.  These 

choices  were  made  because  there  is  a  large  frequency  difference  between 

modes  14  and  15  and  another  large  difference  between  modes  22  and  23. 

However,  it  soon  became  obvious  that  modes  11,  12,  16,  17,  19,  20,  24,  29 
and  30  were  not  affected  by  control.  These  modes  are  local  to  the  antenna 
and  the  counterweight  arms  and  thus  are  not  affected  by  actuators  on  the 
beam.  This  makes  physical  sense,  since  the  antenna  is  isolated  from  the 
beam  by  the  relativly  massive  tip  package.  These  modes  were  put  into  the 
residuals.  It  should  be  noted  that  modes  11  and  12  are  fairly  important  in 
the  motion  of  the  dish  rim  motion,  but  are  nearly  unobservable  and 
uncontrollable  from  anywhere  else  on  the  structure,  including  the  dish 
attachment  point.  Similarly,  modes  19  and  20  only  affect  the  motion  of  the 


couterweight  arms.  Modes  16  and  17  are  less  local,  but  nontheless  are  only 
slightly  affected  by  control.  Modes  29  and  30  are  local  to  the  antenna  dish, 
and  are  not  touched  by  the  controller. 

Also  immediately  obvious  after  a  few  runs  was  that  mode  28  was  very 
important  to  system  stability.  Thus  it  was  made  into  a  suppressed  mode. 
After  finer  adjustments,  mode  27  was  suppressed.  It  was  very  interesting  that 
although  both  modes  27  and  28  became  much  more  stable  when  treated  as 
residual,  other  modes  lost  stability  unless  these  two  modes  were  suppressed. 
Similarly,  the  eigenvalues  of  modes  22,  23,  and  25  moved  significantly  when 
treated  as  residual,  so  these  modes  were  suppressed.  As  has  been  discussed 


previously,  suppression  lowers  the  observability  of  the  controlled  modes.  This 
was  evident  in  three  cases.  Mode  13  is  a  torsion/x  bending  mode,  as  is 
mode  15.  Many  runs  were  made  with  mode  15  suppressed  and  for  these 
runs,  mode  13  lost  damping  and  was  nearly  unstable.  Once  mode  15  was 
moved  into  the  controller,  mode  13  became  controllable.  Similarly,  modes  14 
and  18  are  both  y  bending  modes.  If  mode  18  is  suppressed,  mode  14  gains 
very  little  damping  by  being  controlled.  Thus  mode  18  was  moved  into  the 
controller.  The  final  example  of  loss  of  observability,  and  the  most  difficult 
case  is  the  interaction  between  the  first  bending  mode  and  the  z— axis 
rotational  free  body  mode.  These  modes  are  difficult  to  stabilize 

simultaneously.  The  real  problem  is  that  the  mode  4  is  not  very  observable. 
If  the  free  body  mode  is  made  a  residual,  it  goes  unstable  and  the  torsion 
mode  is  stable  with  fairly  good  damping.  If  the  free  body  mode  is 
suppressed,  the  torsion  mode  is  not  controllable.  If  both  modes  are 
controlled,  the  gains  can  be  set  so  that  one  or  the  other  is  controllable,  but 
not  both.  This  problem  can  be  solved  by  dropping  the  z— axis  torque  motor 
from  the  control  loop  and  making  the  z-axis  free  body  mode  a  residual. 
This  is  not  a  very  palatable  solution,  although  a  stable  configuration.  Mode 
4  can  be  adequately  controlled  by  the  other  actuators  in  the  system  since 
there  is  strong  coupling  between  the  torsion  modes  and  the  x  plane  bending 
modes 

Also  a  problem  with  the  controller  is  that  if  the  sensors  at  the 
attachment  to  the  ceiling  are  included  in  the  model,  the  x  and  y  free  body 
modes  are  unstable.  This  is  because  the  only  actuators  that  can  control 
these  modes  are  those  in  the  base  exitation  table  (BET).  Since  these  two 
sensors  don’t  contribute  information  on  any  bending  modes,  they  can  be  left 
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out  of  the  control  loop.  If  they  are  included  in  the  control  loop,  the  BET 
must  also  be  in  the  loop. 

Unfortunately,  using  the  BET  and  leaving  the  sensors  in  the  loop  doesn’t 
contribute  to  the  solution  of  the  problem  with  the  first  torsion  mode.  This 
problem  is  intensified  by  the  fact  that  modes  23  and  25  are  torsion/bending 
modes.  Thus  if  these  are  suppressed,  the  observability  of  all  the  torsion 
modes  is  severely  reduced.  An  inspection  of  the  $t  observation  matrix  shows 
that  mode  4  only  shows  up  significantly  in  two  sensors.  This  demonstrates 
that  it  is  important  to  have  enough  torsion  sensors  with  this  structure. 

The  final  iteration  with  the  available  configuration  has  modes  23  and  25, 
both  torsion/bending  modes,  in  the  controller.  This  allows  mode  4  to  be 
controlled.  The  insight  that  lead  to  this  assignment  was  that  the  4th  mode 
really  only  has  two  sensors  that  measured  it  well,  thus  when  these  modes 
were  suppressed  they  took  all  of  the  obsevability  away  from  mode  4.  It  still 
requires  the  highest  cost  of  any  mode  in  order  to  gain  significant  damping. 
The  other  problem  with  this  is  that  it  is  not  prudent  to  control  a  mode  of 
that  high  a  number  in  the  model.  Higher  frequency  modes  are  normally 
modelled  less  accurately  than  lower  frequency  modes  when  using  finite 
element  techniques.  Figure  5  is  a  plot  of  the  eigenvalues  of  this  control 
configuration  showing  the  eigenvalues  both  before  and  after  suppression. 
Before  control,  all  eigenvalues  are  on  the  straight  line  associated  with  f  = 
0.01,  the  assumed  passive  damping  ratio.  As  can  be  seen  from  figure  5,  the 
eigenvalues  of  the  suppressed  modes  lie  almost  directly  on  this  line.  The 
importance  of  suppression  to  system  stability  can  be  seen  in  all  eigenvalue 
plots  generated  during  this  study.  In  every  case,  the  pre— suppression  system 
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It  was  then  proposed  that  additional  rate  gyros  be  placed  at  the 
momentum  exchanger  positions.  In  application,  this  could  be  implemented 
by  realigning  the  two  accelerometers  that  are  already  at  these  positions, 
taking  sum  and  difference  calculations  to  produce  both  acceleration  and 
angular  information.  Since  the  bottom  proof- mass  actuator  is  very  close  to 
a  rate  gyro,  an  accelerometer  was  placed  at  the  antenna  attachment  point 
pointed  in  the  circumferencial  direction.  With  these  sensors  in  place  modes 
23  and  25  could  be  suppressed  while  still  adequately  controlling  the  first 
torsion  mode.  Figure  6  is  a  plot  of  the  eigenvalues  of  this  control 
configuration. 

Another  aspect  of  this  study  was  to  explore  the  performance  of  a  two 
controller  system.  Since  both  controlability  and  observability  are  given  up  in 
order  to  produce  decoupled  controllers,  it  should  not  be  expected  that  a  two 
controller  system  would  perform  as  well  as  a  one  controller  system  given 
sufficient  computer  speed.  In  fact  the  application  of  two  controllers  to  this 
system  did  result  in  a  stable  system,  but  the  eigenvalues  of  the  system  were 
moved  very  little.  In  order  to  stabilize  the  system,  mode  21  had  to  be 
controlled.  In  the  other  designs,  mode  21  had  been  treated  as  a  residual. 
However,  since  in  the  two  controller  model,  modes  of  higher  frequency  than 
mode  21  were  being  controlled  rather  than  suppressed,  mode  21  was  being 
subjected  to  control  energy  and  thus  went  unstable.  Figure  7  is  a  plot  of 
the  eigenvalues  of  this  control  design. 

The  final  part  of  the  investigation  explored  the  application  of  an  actuator 
to  the  tip  package.  This  was  done  in  order  to  demonstrate  its  desirability, 
not  to  show  that  the  results  could  actually  be  accomplished.  The  actuator 
was  assumed  to  be  of  zero  weight  thus  not  changing  the  mode  shapes.  The 
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The  time  response  of  the  system  to  various  initial  conditions  was 
considered  next.  The  time  response  for  the  antenna  attachment  point  was 
considered  representative  of  the  overall  system  response  and  so  was  plotted 
for  all  cases.  Figure  9  shows  the  coordinate  system  used  for  this  discussion. 
Figures  10  and  11  are  the  response  of  the  17  sensor  and  15  sensor  systems 
to  a  .68  m/s  velocity  in  the  y  direction.  The  curve  with  the  largest 
magnitude  is  the  y  position  and  the  other  curve  is  the  z  position  in  both 
plots.  The  significant  vibration  in  both  plots  occurs  in  the  y  pendulum 
mode.  The  z  position  plots  show  evidence  of  both  the  pendulum  mode  and 
the  first  bending  mode. 

Figures  12  and  13  are  the  response  to  an  initial  .36  m/s  velocity  in  the  x 
direction.  Figure  12  shows  the  response  of  the  17  sensor  model,  and  the 
significant  vibration  is  in  the  x  pendulum  mode  with  some  evidence  of  the 
first  torsion  mode.  The  15  sensor  system  shows  much  more  vibration  in  the 
torsion  mode  with  some  pendulum  mode. 

Figures  14  and  15  show  the  response  of  the  system  to  a  velocity  in  the  x 
plane.  This  again  shows  mainly  vibration  in  the  y  pendulum  mode.  The  15 
sensor  model  again  shows  a  bit  of  vibration  in  the  torsion  mode. 

Figures  16  and  17  show  the  response  of  the  17  and  15  sensor  models  to 
an  initial  rotational  velocity  about  the  z  axis.  Although  the  15  sensor  model 
shows  somewhat  higher  overshoot,  it  quells  the  system  down  much  faster. 
This  is  to  be  expected  because  there  are  more  controlled  modes  in  this 
model  than  in  the  17  sensor  system.  Since  the  observability  of  the  system  is 
degraded  to  suppress  more  modes  in  the  17  sensor  system,  the  15  sensor 
model  should  be  expected  to  perform  better  in  this  situation.  However,  there 
is  more  uncertainty  that  the  simulated  15  sensor  response  could  actually  be 
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achieved. 


Figure  18  is  the  response  to  the  same  initial  conditions  of  the  system 
model  with  the  extra  assumed  actuator  at  the  tip  package  position.  As  can 
be  seen,  this  effectivly  stops  the  system  from  rotating  in  much  less  time  than 
either  of  the  other  one  controller  model. 

Finally,  figure  19  is  the  response  of  the  two  controller  system.  The  system 
response  is  poor  as  would  be  expected  for  this  sensor  actuator  configuration. 
The  two  controller  model  would  have  also  benefitted  from  an  additional 
actuator,  but  this  case  was  not  studied. 
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VII.  Recommendations  and  Conclusions 


Since  the  object  is  to  actually  run  a  control  law  (in  basic)  the  more 
benign  configurations  should  be  studied.  I  feel  confident  that  the  cruciform 
is  controllable  using  direct  output  feedback,  if  only  because  they  have  been 
controlling  it  using  loop  closing  methods. 

Other  practical  issues  are  the  fact  that  we  have  a  barely  stable  system 
under  computer  control.  Since  we  have  designed  a  continuous  time  controller 
assuming  high  enough  bandwidth,  this  should  be  checked.  I  have  worked  out 
the  math  to  digitize  the  system,  and  then  the  eigenvalues  can  be  evaluated 
using  the  same  IMSL  routines  as  presently  used,  with  stability  ensured  if  the 
eigenvalues  are  inside  the  unit  circle.  The  only  element  of  the  digitized 
system  not  presently  implemented  is  the  convolution  integral.  The  MEXP 
subroutine  can  be  used  to  digitize  the  A  matrix.  Since  the  convolution  has 
to  be  done  to  provide  a  forced  system  response,  this  is  a  logical  and 
necessary  extension. 

The  continuous  time  transformation  techniques  have  exact  parallel 
techniques  in  discrete  time.  The  suppression  of  observation  spillover  is 
exactly  the  same.  The  technique  to  suppress  control  spillover  is  exactly  the 
same,  except  that  the  technique  must  be  performed  on  Bd,  the  digitized 
control  input  matrix.  The  only  real  problem  is  that  the  discrete  time  and 
continuous  time  algebraic  Riccati  equations  are  not  the  same  for  t  greater 
than  zero.  The  LQGLIB  library  available  on  the  ASD  Cyber  has  routines 
that  will  solve  the  discrete  time  Riccati  equation. 

One  final  issue  that  must  be  resolved  is  that  of  real  actuator  dynamics. 
This  and  past  thesis  efforts  have  treated  actuators  as  perfect  devices  with  no 
dynamics.  This  is  probably  fairly  good  for  the  torque  motors,  but  proof 


mass  actators  have  first  order  dynamics  at  best,  and  are  nonlinear  in  reality. 
TRW  has  a  non-linear  actuator  truth  model.  Some  actuator  dynamics  would 
be  included  in  a  reasonable  performance  evaluation.  The  incorporation  of 
actuator  dynamics  into  the  system  model  in  order  to  suppress  and  control 
means  resolving  the  eigenvalue-eigenvector  problem. 


Appendix  A 
Program  Source  Code 
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PROGRAM  ACONE 

C  ATEST1  WITH  ONE  TRANSFORMATION 
C  THIS  PROGRAM  IS  A  CLONE  OF  D0FB33 
C  IT  ONLY  USES  ONE  CONTROLLER 

C  THIS  PROGRAM  IS  MODIFIED  TO  WORK  WITH  THE  NASA  GROUND  TEST  VEHICLE 
C  IT  USES  ONE  CONTROLLER ,  ALTHOUGH  IN  THE  FUTURE  MORE  MAY  BE 
C  AVAILABLE 
C 

C  TO  CHANGE  THE  PROBLEM  TO  OTHER  CONFIGURATIONS 
C  THE  DIMENSION  OF  ALL  THESE  ARRAYS  HAS  TO  BE  GREATER  THAN 
C  MAX (  #  OF  SENSORS,  #  OF  ACTUATORS,  *  OF  MODES  PER  CONTROLLER) 

C  UNFORTUNATELY,  YOU  MUST  ALSO  DIMENSION  MAJM  EXACTLY  THE 
C  SAME  SIZE  AS  YOUR  PROBLEM. 

C  OTHER  unfortunates: 

C  ABCG  DESTROYS  THE  K  MATRIX 
C  FINDK  IS  HARDWIRED  FOR  NCOL  =  41 

C  IF  YOU  HAVE  MORE  THAN  *0  MODES  IN  THE  CONTROLLER,  SUPPRESSED  OF 
C  IN  THE  RESIDUALS  YOU  WILL  HAVE  TO  MODIFY  FINDK 
C 

C  THAT  IS  ALL  BESIDES  THE  INITIALIZATION  SECTION  AFTER 
C  THE  VARIABLE  DECLARATIONS.  (SEE  THE  FILE  ASTART  FOR  A  DISCUSSION 
C  OF  WHY  THESE  CAN  BE  DIMENSIONED  THIS  WAY.) 

C  OTHERWISE,  THE  PROGRAM  GETS  ITS  INFO  FROM  THE  INPUT  FILE 
C 

C  NEXT  THREE  LINES/  SYSTEM  PARTITIONS 

C  DIMENSIONED  NCOL  BY  NCOL 

REAL  A1(41,41),A2(41,41),A3<41,41) 

REAL  B 1(41, 41), B2 <41, 41), B3 (41,41) 

REAL  Cl (41 ,41 ) ,C2(41 ,41 ) ,C3(41 ,41 ) 

C  QA1  »  WEIGHTING  MATRIX  FOR  THE  CONTROLLER 
REAL  0A1 (41,41) 

C  THESE  THREE  MATRICES  ARE  AFTER  SUPPRESSION 
REAL  CSTAR1 (41,41) 

REAL  GT1 (41 ,41 ) 

REAL  KHAT1 (41 ,41) ,KSTAR1(41 ,41) ,KW0RK(41,41 ) 

REAL  SAT (41 ,41 ) ,AKC(41 ,41) , ACT (41 ,41 ), BCG (41 ,41 ) 

REAL  P(41,41)  ,S(41,4.l )  ,KCC(41 ,41 ) 

REAL  ABG1 (41 ,41 ) ,ABG2(41 ,41 ) ,ABG3(41 ,41 ) 

REAL  ABG5 (41,41) , ABG6( 41 ,41 ) ,ABG7(41 ,41 ) 

REAL  GAIN1 (41 ,41 ) ,KT1 (41 ,41 ) 

REAL  GAMMAK41 ,41)  ,T1  (41 ,41 )  ,TKG1  (41 ,41 ) 

REAL  TRT(41 ,41)  ,TEN(41 ,41 )  ,CT<  41 ,41 )  ,V(41,4.l ) 

REAL  MAJM (60,60) , D ( 33 ) , XO ( 60 ) , W( 33 ) , TOL , DT , XC( 60 ) 

REAL  Y ( 41 ) 

REAL  ZETA,AA(33) ,BB(33) ,SING( 41 ) ,XTR( 41 , 41 ) ,X1 ( 60) 

REAL  EAT2(60,60> , WORK (60,60 ) ,STOR( 41, 41) 

REAL  PHI A (41 ,41) , PHIS (4 1,41) , PHILOS (3, 60) ,INIT<2,33) 


55 


non  onnnnn  noonnnnoon 


I NTEGER  N , N2 , NCI , NSUP , NC 12 , NSUP2 , NRES , PUT 
INTEGER  I , J,K,L,M,KK,LL,MM 

INTEGER  CJ,NACT,NSEN,  IRC  41 )  ,  I ER, SNIP, NCOL, NCOL 1 
INTEGER  NDA,NDIM,NDA1 ,NDIM1 ,ZZ,E1 ,P1 
COMPLEX  Z ( 66 ) , W1 ( 33 ) 

C0MM0N/MAINA/NDA,NDA1 

COMMON/M A I NB/NCOL  , NCOL 1 

COMMON/MAINl/NDIM,NDIMl, TEN, X( 3700) 

COMMON/MAIN2/STOR 

C0MM0N/MAIN3/XTR 

COMMON/SAVE/T( 100) ,TS< 100) 

COMMON/ INOUT/KOUT, TAPE 
COMMON/ I O/NSEN, N ACT 
COMMON/PRO/Q 

COMMON/ NUM/IC1 (41),IC2(41),IC3(41), IR,NC1 , NSUP, NRES 


INITIALIZATIONS 

THE  NEXT  6  VARIABLES  SHOULD  NOT  BE  CHANGED  UNLESS 
THE  DIMENSION  OF  YOUR  PROBLEM  IS  SUCH  THAT  THE 
NUMBER  OF  MODES  IS  GREATER  THAN  33  OR  YOU  HAVE  MORE  THAN 
A 1  ACTUATORS  OR  SENSORS.  OR  MORE  THAN  15  MODES  IN  A  CONTROLLER. 

NDIM  IS  GREATER  THAN  THE  NUMBER  OF  MODES 

NDIM  =  41 
NDIM1  =  42 

NCOL  IS  THE  LENGTH  OF  A  COLUMN  IN  THE  WORN  MATRICES  SUCH 
AS  PHIA,  PHIS,  Al,  A2 ,  A3,  XETC. 

THESE  ARRAY  SIZES  ARE  NOT  IMPORTANT  AS  LONG  AS  THEY  ARE  BIGGER  THAN 
YOUR 

ACTUAL  PROBLEM  SIZE 
NCOL  =41 
NC0L1  =  42 
NDA  IS  TWICE  NDIM 

AND  THUS  GREATER  THAN  THE  DIMENSION  OF  THE  FIRST  ORDER  SYSTEM 

NDA  =  60 
NDA 1  =  61 

C  DON'T  REALLY  KNOW  WHAT  THESE  DO... 


on  nonnnnnn  no 


OPEN  <  8 , FILE  = ' TAPE8 ' ) 

REWIND  8 

OPEN < 6 , F ILE= ' TAPE6 ' ) 

REWIND  6 
CALL  HEADER  ( ) 

HEADER  PRINTS  THE  HEADER  AT  THE  TOP  OF  THE  OUTPUT 
XZZZXZZZZZZZZZZZZZZXZXZZZZZZZZZZXXZZZXZZZXZZZZZZXX 
READ  IN  INFORMATION  FROM  INPUT  FILE 

zzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzzz 

IF  Q  EQ  2,  THEN  PRNT,  PRNTXL,  PRINTK  DON'T  WORK 
PRINT*,'  TO  PRINT  ALL  OF  THE  MATRICES  ENTER  1,  ELSE  ENTER  0 
READ<8,*)  Q 
PRINT*, Q 

PHI  MATRICES  AND  CONTROLLER  ENTRIES 
PRINT' (///) ' 

NCI  IS  THE  NUMBER  OF  CONTROLLED  MODES 
NSUP  IS  THE  NUMBER  OF  SUPPRESSED  MODES 
NRES  IS  THE  NUMBER  OF  RESIDUAL  MODES 
NACT  IS  THE  NUMBER  OF  ACTUATORS 
NSEN  IS  THE  NUMBER  OF  SENSORS 

ZETA  IS  A  COMMON  ASSUMED  DAMPING  RATIO  FOR  ALL  MODES 

PRINT*,'  ENTER  NCI , NSUP, NRES, NACT, NSEN, ZETA  >' 

READ  <  8 , * )  NCI, NSUP , NRES , NACT , NSEN , ZETA 
PRINT*, NCI , NSUP, NRES, NACT, NSEN, ZETA 
PRINT*,'  ' 

PRINT*,'  ENTER  THE  ' , NACT , '  ELEMENTS  FOR  EACH  PHIA  ' 
PRINT*,'  ' 

N  =  NCI  +  NSUP  +  NRES 


PRINT*,'  ENTER  PHIA,  A  ',N,'  BY  ' , NACT , '  MATRIX 
*  BY  ROWS' 

DO  1  1=1, N 

READ (8,*)  (PHIA(I,J) ,J=1,NACT> 

CONTINUE 

CALL  PRNT ( PH I A, N, NACT) 

PRINT' <//) ' 

PRINT*,'  ENTER  THE  ' , NSEN , '  ELEMENTS  FOR  EACH  PHIS  ' 
PRINT*,'  ' 

PRINT*,'  ENTER  PHIS,  A  '.NSEN,'  BY  ',N.'  MATRIX 


on  o  n  o 


*  BY  COLUMNS' 

C  THE  PHIS  MATRIX  IS  FROM  THE  EQUATION 
C  Y  =  C  *  X  WHERE  C  IS  GIVEN  BY 
C  C  =  <  0  :  PHIS  TRANSPOSED  ) 

C  IN  THIS  PROGRAM,  BECAUSE  ONLY  VELOCITY 
C  SENSORS  ARE  USED .  THIS  IS  A  FUNCTION 
C  OF  SUBROUTINE  FORMC 
DO  2  1=1, N 

READ ( 8 , *  >  (PHIStI, J) , J=1,NSEN) 

2  CONTINUE 

CALL  PRNT  <PHIS,N,NSEN) 

PRINT' (//) ' 

OMEGAS 

PRINT*,  '  ENTER  THE  VALUE  FOR  EACH  OMEGA  ' 
PRINT*,'  ' 

READ (8,*) <W< I ) , 1=1 ,N> 

DO  3  1  =  1,  N 

PRINT*,'  ' , W ( I ) 

D < I )  =  -2.  *  ZETA  *  U(I) 

3  CONTINUE 


20  CONTINUE 
C 

PRINT' <//)  ' 

C 

C  MODE  ASSIGNMENT  TO  CONTROLLERS 

C 

c 

PRINT' (///)  ' 

PRINT*,'  THE  FOLLOWING  MODES  ARE  ENTERED  ACCORDING  TO  THE 
PRINT*,'  ORDER  IN  WHICH  THEY  ARE  ENTERED  IN  THE  DATA  FILE 
PRINT*,'  AND  NOT  ACCORDING  TO  THEIR  ACTUAL  MODE  NUMBER.  ' 
PRINT' <//> ' 

PRINT*,'  ENTER  THE  ',NC1,'  CONTROLLED  M(  :S  >' 

READ  <  8 , * )  <IC1(I),I  =  1, NC 1 ) 

PRINT*,'  ' , (IC1<I) ,1=1, NCI) 

PRINT*,'  ' 

PRINT*,'  ENTER  THE  '  ,NSUF’,  '  SUPPRESSED  MODES  >' 

READ  <  8 , *  >  < IC2< I > , 1  =  1 , NSUP ) 

PRINT*,'  ', (IC2(I) ,1=1, NSUP) 

PRINT*,'  ' 

PRINT*,'  ENTER  THE  ' , NRES , '  RESIDUAL  MODES  >' 

READ (3,*)  ( IC3 ( I > , I =1 , NRES ) 

PRINT*,'  ' , (IC3( I) ,1=1 ,NRES) 

PRINT*,'  ' 
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PRINT*,'  ' 


NC12  =  2  *  NCI 
NSUP2  =  2  *  NSUP 
NRES2  =  2  *  NRES 
N2  -  2  *  N 

M  =  NC12  +  NSUP2  +  NRES2 
NDA  =  M 
NDA1  =  M  +  1 
NDIM  =  N 
NDIM1  =  N  +  1 
100  CONTINUE 

PRINT' (///)' 

READ  IN  THE  WEIGHTING  MATRIX  DIAGONAL 
VALUE  FOR  EACH  MODE 

PRINT*,'  ENTER  THE  DIAGONAL  VALUES,  IN  MODE  INPUT 
PRINT*,'  ORDER,  FOR  THE  CONTROL  WEIGHTING  MATRIX 
READ ( 8 , * )  (AA(IQ.l)  THEN 

PRINT*,'  THE  CONTROLLER  A  MATRIX  tS  ' 

CALL  PRNT ( A1 ,NC12 ,NC12) 

PRINT*,'  THE  CONTROLLER  B  MATRIX  (S  ' 

CALL  PRNT  <  B1 ,NC12 , NACT ) 

PRINT*,'  THE  CONTROLLER  C  MATRIX  IS  ' 

CALL  PRNT (Cl ,NSEN,NC12> 

PRINT*,'  THE  Cl  CONTROL  WEIGHTING  MATRIX  IS  ' 
CALL  PRNT ( QA1 ,NC12 » NC12) 

PRINT*,'  THE  SUPPRSD  MODES  A  MATRIX  IS  ' 

CALL  PRNT  ( A2 ,  NSUP2 ,  NSUP2  ) 

PRINT*,'  THE  SUPPRSD  MODES  B  MATRIX  IS  ' 

CALL  PRNT <B2,NSUP2, NACT) 

PRINT*,'  THE  SUPPRSD  MODES  C  MATRIX  IS  ' 

CALL  PRNT ( C2 , NSEN , NSUP2 ) 

IF  (NRES.GT.O)  THEN 

PRINT*,'  THE  RESIDUAL  A  MATRIX  IS  ' 

CALL  PRNT ( A3 , NRES2 , NRES2 ) 

PRINT*,'  THE  RESIDUAL  B  MATRIX  IS  ' 

CALL  PRNT ( B3 ,NRES2 , NACT ) 

PRINT*,'  THE  RESIDUAL  C  MATRIX  IS  ' 

CALL  PRNT (C3, NSEN, NRES2) 

FI.SE 

PRINT*,'  NO  RESIDUAL  TERMS  ' 

END  IF 
END  IF 
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C  THIS  SECTION  GENERATES  THE  R.ICCATI  SOLUTION 
C  AND  THE  GAIN  MATRIX  FOR  THE  CONTROLLER 

C 

CALL  VMULFP <  B1 , B1 , NC12 , NACT , NC 12 , NCOL , NCOL , SAT , NCOL ,  IER ) 
IER  =  0 
TOL  =  O.OOl 

PRINT*,'  THE  FOLLOWING  ARE  THE  MRIC  A+BG  1  INPUTS  ' 
PRINT' (//) ' 

PRINT*,'  THE  MATRIX  A1  IS  ' 

CALL  PRNT (A1,NC12,NC12) 

PRINT*,'  THE  MATRIX  SAT  <B1*B1T)  IS  ' 

CALL  PRNT ( SAT , NC12 , NC12 ) 

PRINT*,'  THE  MATRIX  QA1  IS  ' 

CALL  PRNT ( QA1 , NC12 , NC12 ) 

PRINT*, '  NC12  =  '  , NC12 
PRINT' <//) ' 

CALL  MRIC < NC12, A1 , SAT, QA1 ,S,ABG1 , TOL, IER) 


C 

c 

c 


629 


ABG1  =  A.l  +  B1G1 


PRINT*,'  THE  EIGENVALUES  OF  A1  +  B1G1  ' 

PRINT' </) ' 

CALL  EIGRF  <  ABG1 , NC12 , NCOL ,0»W1,TEN, NCOL , STOR , IER  > 
DO  629  I  =  .1 ,  NC  1 2 
PRINT*,'  '  , U1 < I ) 

PRINT' </)  ' 


IN  A  MULTIPLE  CONTROLLER,  SOME  OF  THE  RICCATI  GAINS 
HAVE  TO  BE  RECOMPUTED  USING  THE  TRANSFORMATION  MATRICES 


C 
C 
C 

C  THIS  IS  NOT  SO  FOR  ONE  CONTROLLER  BECAUSE  WE  SUPPRESS 
C 

c 
c 


THE 

THE 


C  MATRIX,  AND  THE  C  MATRIX  IS  NOT  INVOLVED  IN 
RICCATI  SOLUTION. 


C 

C 

C 

C 

C 

C 

C 


PRINT*,'  THE  RICCATI  SOLUTION  OF  AC  +  BCG  *1  IS  ' 

CALL  PRNT (S,NC12,NC12) 

IN  ORDER  TO  HAVE  A  CONTROL  WIEGHTING  MATRIX  R,  THIS  CODE  WOULD  BE 


CHANGED.  G  =  -RA-1*BAT*S,  BUT  OUR  R 


I  UNTIL  CODED  DIFFERENTLY 


CALL  VMULFM <  B1 , S , NCI 2 , NACT , NCI 2 , NCOL , NCOL , GAIN1 , NCOL , IER ) 
PRINT*,'  THE  G1  GAIN  MATRIX  IS  ' 

CALL  PRNT ( GAIN1 , NACT, NC12) 

THIS  IF  STATEMENT  GETS  RID  OF  THE  UNSUPPRESSED  EIGENVALUE 
CALCULATIONS.  AFTER  THE  PROGRAM  WORKS,  THESE  REALLY  LOSE 
MOST  OF  THEIR  SIGNIFICANCE,  SO  WHY  WASTE  THE  TIME? 

IF  (Q.LT.3)  THEN 

FINDK  FINDS  THE  K  GAIN  MATRIX  FROM  THE  OPIMAL  FULL  STATE 


ul*-'.  JU.-J 

.*•  . 


.  ,4 


- 


*  O-  . 
'  *■*  .  ' 
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non  o  o  ononnnnnonnnonoonoo  o  on 


FEEDBACK  GAIN,  G. 

CALL  FINDK ( Cl , NSEN , NC12 , WORK , GAIN1 , NACT , KHAT 1 ) 
PRINT#, ' THE  FEEDBACK  MATRIX  KHAT 1  IS' 

CALL  PRNT ( KHAT 1 , NACT, NSEN) 

DO  13  1=1, NACT 
DO  13  J=1 ,NSEN 
13  KWORK ( I , J )  =  KHAT 1 < I , J ) 


THIS  SECTION  GENERATES  THE  BLOCK  SEGMENTS 
OF  MAJM  AND  PUTS  THEM  INTO  THE  MAJM  MATRIX 


THE  SYSTEM  MATRIX  LOOKS  LIKE: 


**** 

***# 

* 

* 

*  A1+B1K1C1 

B1K1C2 

B1K1C3 

* 

* 

* 

*  B2KIC1 

A2+B2KIC2 

B2KIC3 

* 

* 

* 

*  B3KIC1 

B3KIC2 

A3+B3KIC3 

* 

* 

* 

**** 

*#** 

K=NC12 
L-=NSUP2  +  K 

MM  =  NCI  2  +  NSUF‘2  +  NRES2 

PRINT ' <///>' 

DO  16  1=1, NACT 
DO  16  J=1 , NSEN 

16  KHAT 1 ( I , J)  =  KWORK< I , J) 

DO  200  1=1, MM 
DO  200  J= 1 , MM 
200  MAJM ( I , J )  =0.0 

NOTE  THAT  ABGC  DESTROYS  KHAT 1  EACH  TIME!!! 

CALL  ABGC ( A1 ,B1 , Cl , NCI 2, KHAT 1 , NACT, NSEN, ABG1 ) 

DO  17  1=1, NACT 
DO  17  J=1 ,NSEN 

17  KHAT 1 ( I , J )  =  KWORK ( I, J) 


noon  on  non  ooo  noon  no 


PARTITION  1,1 

HO  201  1=1 , NC12 
DO  201  J=1 ,NC12 

201  MA  JM ( I , J )  =  ABG1 ( I , J ) 

PARTITION  2,2 

CALL  ABGC <  A2 , B2 , C2 , NSUP2 , NHAT 1 , NACT , NSEN , ABG2 ) 

HO  18  1=1, NACT 
HO  18  J=1 , NSEN 
18  NHAT 1 < I , J )  =  KWORK  < I , J ) 

HO  202  1=1 ,NSUP2 
no  202  J=1 ,NSUP2 

202  MAJM< I+K, J+K)  =  ABG2(I,J) 

PARTITION  3,3 

CALL  ABGC <  A3 , B3 , C3 , NRES2 , NHAT 1 , NACT , NSEN , ABG3 ) 
HO  24  1=1, NACT 
DO  24  J=1 , NSEN 
24  NHAT 1 < I , J )  =  NWORN ( I , J) 

HO  203  1=1 ,NRES2 
HO  203  J=1,NRES2 

203  MAJM<  I+L,  J+l.)  =  ABG3(I,J) 

PARTITION  1,2 

CALL  MMUL (Bl, NHAT 1 ,NC12, NACT, NSEN, TEN) 

CALL  MMUL ( TEN , C2 , NC 1 2 , NSEN , NSUP2 , BCG ) 

HO  208  1=1 , NC12 
HO  208  J=1 ,NSUP2 
208  MAJM ( I , J+N )  =  BCG< I , J ) 

, NRES2 

405  PRINT*,'  '  ,  W1  <  I  > 

PRINT' <//>' 

ELSE 

PRINT*,'  NO  RESIHUAL  TERM  EIGENVALUES  ' 

PRINT' (/) ' 

END  IF 
ENHIF 


THIS  SECTION  FORMS  THE  TRANSFORMATION  MATRIX. 
SINCE  THE  OBJECT  OF  THIS  PROGRAM  IS  TO  BESIGN  ONE 


C  CONTROLLER 

C  U  =  K  *  GAMMA  *  C  *  X 

C 


c 

AFTER 

THE  TRANSFORMATION 

IS  COMPLETE 

, 

c 

THE  MAJM  WILL  LOOK  LIKE: 

c 

c 

***# 

**** 

c 

* 

* 

c 

* 

( A1+B1K1C1 ) 

0 

B1K1CR 

* 

c 

* 

* 

c 

* 

B2K1C1 

'  A2 

B2K1CR 

* 

c 

* 

* 

c 

* 

(BRK1C1 ) 

( BRK1C2 ) 

( AR+BRK1CR ) 

* 

c 

* 

* 

c 

**** 

**** 

c 

c 

C  WHERE  THE  NON-ZERO  TERMS  INCLUDE  THE 
C  TRANSFORMATION  MATRIX. 

C 

C  GENERATE  THE  TRANSFORMATION  MATRIX  GAMMA1 
C 

CALL  TFR ( CT ,  C2 , NSEN , NSUP2 ,1,2) 

DO  500  1=1, NSUP 
DO  500  J=1 , NSEN 
500  0(1, J)  =  CT  < I +NSUF , J ) 

NRV  =  NSUP 
PRINT*,'  V  ( C2)  IS  ' 

CALL  F'RNT(V, NRV, NSEN) 

CALL  LSVDF ( V , NCOL , NRV , NSEN , TEN , NCOL , -1 , SING , STOR , I ER  > 
PRINT*,'  ' 

PRINT*,'  V  OUT  OF  LSVDF  IS  ' 

CALL  PRNT(V, NSEN, NSEN) 

PI  =  NSEN  -  NRV 
IF  (Pl.LT.l)  THEN 
DO  503  1=1, NSEN 

503  GAMMA1 (1,1)  =  V(I,NSEN) 

PI  =  1 

ELSE 

DO  504  1=1, NSEN 
DO  504  J=1,P1 

504  GAMMA1 ( I , J )  =  V(I,J+NRV) 

END  IF 

C 

PRINT*,'  TRANSFORMATION  MATRIX  GAMMA1  ' 

CALL  PRNT ( GAMMA1 , NSEN, PI ) 

C 

C  CHECK  TO  SEE  THAT  GAMMA1  IS  ORTHOGONAL  TO  C2 
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note:  akc  in  this  section  is  just  a  work  area  to  test 

THE  ORTHOGONALITY  OF  CT  #  TR.  IN  ALL  CASES  IT 
SHOULD  BE  A  BLOCK  ZERO  MATRIX. 

CALL  TFR(CT,C2,NSEN,NSUP2, 1 ,2) 

CALL  MMUL <CT, GAMMA  1  ,NSUP2,NSEN,F*1  ,AKC> 

PRINT#, '  C2T  #  GAMMA1  ' 

CALL  PRNT (AKC,NSUP2,P1 ) 

PRINT#,'  C2  SINGULAR  VALUES  ' 

CALL  PRNT ( SING , NRV , 1 ) 

DO  411  1=1 , NSEN 
DO  411  J=1,P1 
411  GT1 <  J , I ) =GAMMA1 < I , J ) 

PRINT#,'  GAMMA1  TRANSPOSE  ' 

CALL  PRNT ( GT1 , FI , NSEN ) 

CSTAR  =  GAMMA  #  C 

CALL  MMUL(GT1 ,C1 ,P1 ,NSEN,NC12,CSTAR1 ) 

PRINT#,'  (GT1XC1)  ' 

CALL  PRNT ( CSTAR 1 ,P1 ,NC12) 

GSTAR1  =  GAIN1 

CALCULATE  UT,  V,  QPLUS  WHERE  CSTAR  =  U*Q#VT  AND 

CSTAR+  =  V*GPLUS*UT 
KSTAR  =  GSTAR*CSTAR+ 

CALL  FINDK ( CSTAR1 , P 1 , NC12 , WORK , GAIN1 , NACT , KSTAR 1 ) 
PRINT#,'  KSTAR1 
CALL  PRNT  <  KSTAR 1 , NACT , PI ) 


FORM  TKG  =  KSTAR  #  GAMMA 

CALL  MMUL ( KSTAR1 ,GT1 , NACT, PI , NSEN, TKG1 ) 

PRINT#,'  (KSTAR1 ) (GT1 ) 

CALL  PRNT ( TKG1 , NACT, NSEN) 

C 

DO  348  1=1, NACT 
DO  348  J=1,NSEN 
348  SAT ( I , J ) =0 , 0 
C 

C  KT1  =  -TKG 

C  REMEMBER  THAT  ABGC  TAKES  CARE  OF  THIS  MINUS  SIGN  FOR 
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C  THE  DIAGONAL  BLOCKS 
C 

CALL  UB0AT<SAT,TKG1 , NACT , NSEN , KT 1 ) 


REMEMBER  THAT  ABGC  DESTROYS  TKG1 
DO  2.1  1  =  1,  NACT 
DO  21  J=1 ,  NSEN 
21  KWORK ( I »  J )  =  TKG1 <  I ,  J ) 

DO  350  1  =  1,  M 
DO  350  J  =  1,M 

350  MAJM< I , J>=0.0 

FILL  IN  PARTITION  1 , 1  OF  MAJM 

CALL  ABGC <A1,B1, Cl, NC12, TKG1 , NACT, NSEN, ABG5) 
PRINT*,'  A1  +  BKTKGDCl  ' 

CALL  PRNT ( ABG5 , NC12 ,NC12 ) 

DO  351  1=1 ,NC12 
DO  351  J=1 , NC12 

351  MAJM< I,J)=ABG5(I,J) 

PART  IT  I  '•  2 


DO  22  1=1, NACT 
DO  22  J= 1 , NSEN 
22  TKGl(IfJ)  =  KUORK  < I , J  > 

CALL  ABGC ( A2 , B2 , C2 , NSUP2 , TKG1 , NACT , NSEN , AD' '6 ) 
PRINT*,'  A2  +  B2 <  TKGl )C2  ' 

CALL  PRNT ( ABG6 , NSUP2 , NSUP2 ) 

J=1 , NRES2 

358  MAJM( I , J+L)=KCC( I  ,  J) 

END  IF 

((<((<(<((<<(((<((((((((((((((<((((((((<< 

NEXT,  FORM  E  TO  THE  AT 
BECAUSE  EIGRF  DESTROYS  MAJM 

IF  (SKIP, EG. 1 )  THEN 
TOL  *  0.001 

CALL  MEXP  <  MM , MAJM , DT , EAT2 ) 

PRINT*,'  MEXPOUT  ' 

PRINT' (//> ' 

PRINT*,'  THE  SOLUTION  EAT2  IS  ' 

CALL  PRNTXL.  ( EAT2 ,  MM ,  MM  ) 

END  IF 


C  EAT2  IS  MOW  THE  SOLUTION  10  E  TO  THE  AT 

C 

c 

C  ))))))>)))))))))))))))))))))))))))))))))))))))) 

c 

PRINT*,'  LOWER  TRIANGLE  ONE  CONTROLLER  MAJOR  MATRIX 

PRINT* , '  AFTER  SUPPRESSION,  WITH  RESIDUALS' 

CALL  PRNT  <  MA JM , M , M ) 

PRINT*, ' (///)' 

PRINT*,'  TRANSFORMED  SYSTEM  EIGENVALUES  ' 

PRINT*,'  AFTER  SUPPRESSION' 

CALL  E I GRF <  MA JM , M , NDA , 0 , Z , TEN , NCOL , WORK , I ER ) 

PRINT' </) ' 

DO  366  I  =  1 ,  M 
366  PRINT*,'  ' , Z ( I ) 

PRINT' <//) ' 

PRINT*,'  EIGENVALUES  OF  A1  +  ( B1 ) ( TKG1 ) ( Cl )  ' 

CALL  E IGRF  < ABG5 , NC12 , NCOL , 0 , W1 , TEN , NCOL , STOR , IER ) 
PRINT'!/) ' 

DO  723  1=1 ,NC12 

723  PRINT*,'  ' , W1 < I ) 

PRINT' (/)  ' 

PRINT*,'  EIGENVALUES  OF  A2  +  ( B2 ) ( TKG1 ) ( C2 )  ' 

CALL  El  GRF  <  ABG6 ,  NSUP2 ,  NCOL  ,  0 ,  W1. ,  TEN ,  NCOL  ,  STOR ,  I  ER ) 
PRINT' </)  ' 

DO  724  1=1 , NSUP2 

724  PRINT*,'  ' , W1 ( I ) 

PRINT' </> ' 

C 

IF  (NRES.GT.O)  THEN 

PRINT*,'  EIGENVALUES  OF  A3  +  < B3 ) ( TKG1 ) ( C3 ) 

CALL  E I  GRF  (ABG7,NRES2,  NCOL,  0,  W 1. ,  TEN,  NCOL,  STOR,  IER) 
PRINT' (/)  ' 

DO  725  1=1 , NRES2 

725  PRINT*,'  ' , W1 < I ) 

END  IF 

C 

IF  ( SKIP »NE . 1 )  THEN 

PRINT*,'  NO  TIME  RESPONSE  PRINTOUT  ' 

ELSE 

PRINT' (///) ' 

PRINT*,'  STATE  TRANSITION  MATRIX >(2' 

CALL  PRNT ( EAT2 , MM , MM ) 

CALL  TIMEX < EAT2 , MM , D T , XO , PDT , TMAX , X 1 , XC , PHILOS , 

*IC1 , IC2, IC3,NC1 ,NSUP,NRES) 

ENDIF 


STM  IS  NOW  THE  SOLUTION  TO  ZDOT  =  MAJM  *  Z 
WE  NOW  PROPAGATE  THE  STATE  IN  DT  STEPS 
END  OF  TIME  RESPONSE  SECTION 

PRINT*, '  ************************************* 

PRINT*,'  **  ** 

PRINT*,'  **  END  D0FB335  PROGRAM  ** 

PRINT*,'  **  ** 

PRINT*, '  ************************************* 

PRINT' <///) ' 

END 


START  OF  SUBROUTINES 

))))))))))))))))))))))))))))))))))))))))))))))) 

ABGC 

ABGC  IS  A  BAD  CITIZEN,  IT  DESTROYS  KCC,  TAKE  HEED 

))))))))))))))))))))))))))))))))>)))))))))))))) 

SUBROUT I NE  ABGC  <  A , B , C , N2 , KCC , NACT , NSEN , ABG ) 
COMMON/MAINB/NCOL , NCOL I 

REAL  B ( NCOL , 1 ) , A ( NCOL , 1 ) , C ( NCOL , 1 ) , KCC ( NCOL , 1 ) , ABG ( NCOL . 1 > 

INTEGER  N2, NACT, NSEN 

CALL  MMUL • B , KCC , N2 , NACT , NSEN , ABG ) 

CALL  MMUL< ABG, C,N2, NSEN, N2, KCC) 

CALL  UBOAT (A,KCC,N2,N2, ABG ) 

RETURN 

END 

))))))))))))))))))))))))))))))))))))))))))))))) 

FUNCTION  DOT 

))))))))))))))))))>)))))))))>)>)))>)))))))))))) 


FUNCTION  DOT (NR, A, B) 
DIMENSION  A< 1 ) ,B< 1 ) 
D0T=0, 

DO  1  1=1, NR 

DOT  =  DOT+A< I )*B< I ) 

RETURN 

END 


onnoon  noono 


FINDS  N  FROM  THE  RICCATI  SOLUTION ,  G 
))))>))))))))))))))))>)))))>>>)>))))>)))))))))) 


SUBROUTINE  FINDK(C,M,N,W,G,L,K> 

M=NSEN  N=NC12  W=WORK  SPACE  L=NACT  (G  IS  NACT  BY  NC12) 
COMMON/MA I NA/NDA , NDA 1 
COMMON/MAINB/NCOL ,NCOLl 
REAL  UT<41,41>,VP<41,41),QPLUS<41,41> 

REAL  SI<42) ,CPLUS<41,41> , OUT <41, 41) 

REAL  U<41,41) ,C<NCOL,NCQL) ,G<NCOL,NCOL> 

REAL  K  <  NCOL  , NCOL  >  *  W  <  NDA  , NDA ) 

DO  1  1=1, M 
DO  1  J=1 ,  M 

1  UT ( I ,  J)=0*0 
DO  2  1=1, M 

2  UT ( I , I ) =1 « 0 
DO  3  1=1, M 
DO  3  J=1,N 

3  VP<I,J)=C(I,J) 

CALL  LSVDF <UP,NCOL,M»N,UT , NCOL, M, SI ,U,IER) 

DO  4  1=1, N 
DO  4  J=1,M 

4  QPLUS ( I , J ) =0 . 0 
DO  5  1=1, N 

5  IF<SI(I) .GT. 0.00001)  QPLUS ( I,I)=1/SI(I) 

DO  6  1=1, N 

DO  6  J=1 , N 

6  V<I,J)=VP<I,J) 

CALL  MMUL < QPLUS , UT ,N,M,M,QUT) 

CALL  MMUL ( U , QUT , N , N , M , CPLUS ) 

CALL  MMUL<G, CPLUS, L,N,M,K> 

RETURN 

END 

)))))))))))))))))))))))))))))))))))))))))))))))))))) 

F0RMX2 

<<(<<(<<<<<<(<<<<<<<<<<<<<<<<<((<<<<<<<(<<<<<<<<<<< 

SUBROUTINE  F0RMX2<X0, INIT) 

C0MM0N/NUM/IC1 <41),IC2<41),IC3(41),IR(41),NC1 ,NSUP , NRES ,NR 
REAL  X0<41),INIT(2,41) 

INTEGER  M,I,J,K,L 
DO  1  1=1, NCI 
M  =  IC1 ( I ) 

XO(I)  =  INIT(1,M) 

1  X0< I+NC1 )  =  INIT<2,M) 

J  =  NC1*2 
DO  2  1=1 , NSUP 
M  =  IC2  < I ) 

XO(I+J)  = 


INIT  < 1 ,M) 


n  o  n  n  n  o  n 


2  X0(  I  +  J+NSUF’)  =  INIT  (  2 ,  M ) 

K  =  J  +  NSUF**2 
DO  3  I=1,NRES 
M  =  IC3 ( I ) 

XO ( I +K )  -  INIT(1,M) 

3  XO( I+K+NRES)  =  INIT(2,M) 
RETURN 
END 


SUBROUTINE  FACTOR < N , A , S ,MR > 

C  A=S'S 

DIMENSION  A(l), 3(1) 

COMMON/M A I NB/  NCOL,NCOLl 
COMMON/ I NOUT/KOUT 
TOL=l ♦ E-6 
MR=0 

NN=N*NCOL 
TOL 1=0* 

DO  1  1=1 , NN , NC0L1 
R=ABS ( A ( I ) ) 

1  IF  (R.GT.T0L1)  T0L1=R 

T0L1=T0L1*1 .E-12 
11  =  1 

DO  50  1=1, N 
IM1=I-1 

DO  5  JJ=I , NN , NCOL 
5  S<JJ)=0. 

ID=I I+IM1 

R=A< ID) -DOT ( IM1 ,S< II ) ,S( II )  ) 

IF  <ABS(R) *LT . (T0L*A(ID)+T0L1) )  GO  TO  50 
IF  (R)  15,50,20 
15  MR=- 1 

URITE(KOUT, 1000) 

1000  FORMAT ( 37H0TRIED  TO  FACTOR  AN  INDEFINITE  MATRIX 
RETURN 

20  S( ID)=SQRT(R) 

MR=MR+1 

IF  (I.EQ.N)  RETURN 
L "II+NCOL 

DO  25  J J=L , NN , NCOL 
I J=JJ+IM1 

25  S(IJ)=(A(IJ) -DOT (IM1,S(II),S(JJ) ) )/S( ID) 


»  j.«  .  ■  !!.■  r-  i;  i. 
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50  I I = I I +NC0L 

RETURN 
END 
C 

C  >>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 
C 

C  SUBROUTINE  FORMA 

C 

C  )))))))))))))))))))>)))))))))))))))>>)))) 
C 

SUBROUTINE  FORMA < A , D , U , N, N2 , IC > 
COMMON/MAINB/NCOL 

REAL  A ( NCOL , NCOL ) , U  <  NCOL ) , D  < NCOL ) 
INTEGER  IC(NCOL) ,I,J,N,M 
DO  1  1=1, N2 
DO  1  J=1 ,N2 
A ( I , J ) =0 . 0 

1  CONTINUE 
DO  2  1=1, N 
M=  IC(I) 

A( < I+N) , ( I+N) )=D<M) 

A(I»<I+N>>  =  1,0 
A ( ( I +N ) , I )  =  -<W(M)**2> 

2  CONTINUE 
RETURN 


END 

C 


SUBROUTINE  FORMS < B , PHI , N,N2 , NACT , IC ) 
COMMON/MAINB/NCOL 
REAL  B(NCOL, NCOL) , PHI (NCOL, NCOL) 
INTEGER  IC ( N ) , NACT , N , M , I , J 
DO  1  1=1, N2 
DO  1  J=1 ,NACT 
B ( I ,  J )  =  0.0 

1  CONTINUE 
DO  2  1=1, N 
M  =  IC(I) 

DO  2  J=1 , NACT 
B ( ( N+I ) , J )  =  PHI(M, J) 

2  CONTINUE 
RETURN 
END 


70 


V-V-V-‘«'-  *  -  ,‘ti  U  i~‘  1  -  ’  *  *  «J  *■  'Mi 


non  nononnn 


C 

C  )>>)>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>>> 

C 

C  SUBROUTINE  FORMC 

C 

C  )))))))))))))))>>>>)>)>>)>>>))>))))>)))) 

c 

SUBROUT I NE  FORMC ( C , PHIS , N , N2 , NSEN ,IC> 
COMMON/MAINB/NCOL 

REAL  C < NCOL, NCOL ), PHIS (NCOL, NCOL) 

INTEGER  ICCN) ,M,NSEN,N,N2, I , J 
DO  1  1=1, NSEN 
DO  1  J=1 ,N2 
C<I,J)  =  0,0 
1  CONTINUE 

DO  2  1=1, NSEN 
DO  2  J=1,N 
M  =  IC<J) 

C  ( I , NT  J )  =F'HIS <  M ,  I ) 

CONTINUE 
RETURN 
END 

)))))>)))))))))))))))))))))))))))))>>>)>))))>> 


FORMO 


)))))))))))))))))))))))))))))))))>)))))>)>)))> 

SUBROUTINE  FORMCH  0  ,  A  ,  N ,  IC  > 

COMMON/MAINB/NCOL 

REAL  A(NCOL) , 0 < NCOL , NCOL ) 

INTEGER  I » J,K,M,N,N2, IC(NCOL) 

N2  =  N  #  2 
DO  1  1=1, N2 
DO  1  J=1,N2 
Q ( I ,  J )  =  0,0 

1  CONTINUE 
DO  2  1  =  1,  N 
M  =  IC(I) 

CHI, I)  =  ACM) 

CHI+N,  I+N)  =  Q  < I , I > 

2  CONTINUE 
RETURN 
END 


))))))))))) 


)  )  ) 
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SUBROUTINE  GMINU 


>>)>)))>)>))>)))))>)>)))))))))))))))))>)) 

3UBR0UT I NE  GM I NO ( NR , NC , A , U , MR , HT ) 
DIMENSION  A  <  1 )  ,  U  <  1 ) 

COMMON/MAIN1/  NDIM , NDIM1 , S ( 1 ) 
COMMQN/MAINB/NCOL  , NCOLI 
COMMON/ I NOUT/KOUT 
TOL=l .  E-12 
MR=NC 
NRM1=NR-1 
TOLl=l . E-20 
JJ=1 

DO  100  J=1,NC 
F  AC=DOT ( NR , A  <  J J  > ,A< JJ>> 

JM1=J-1 
JRM=JJ+NRM1 
JCM=JJ+JM1 
DO  20  I=JJ,JCM 
20  U<I)=0. 

U< JCM>=1 .0 

IF  (J.EQ.l)  GO  TO  54 
KK=1 

DO  30  K=1 ,  JM1 

IF  <S(K> .EQ.1.0)  GO  TO  30 

TEMP=-DOT<NR,A< JJ) ,A(KK) ) 

CALL  VADD ( K , TEMP , U  <  J J ) , U  <  KK ) ) 

30  KN=KK+NCOL 
DO  50  L  = 1 ,2 
KK  =  1 

DO  50  K=1 , JM1 
IF  <  S  <  K ) . EQ ♦  0  •  )  GO  TO  50 
TEMP=-DOT  <  NR , A  <  J J ) , A ( KK  >  > 

CALL  0ADD<NR,TEMP,A< JJ) , A  <  KK  >  > 

CALL  VADD<K,TEMP,U< JJ) ,U<KK)  > 

50  KK=KK+NDA 

T0L1=T0L*FAC 

FAC  =  DOT <  NR , A  <  J J  > , A ( J J  >  ) 

54  IF  <FAC.GT.TOLD  GO  TO  70 
DO  55  I =J J , JRM 

55  A  < I ) =0 . 

S< J)=0. 

KK=  1 

DO  65  K=1,JM1 

IF  ( S < K ) . EQ . 0 , )  GO  TO  65 

TEMP=-DOT<K,lKKK>  ,U<  JJ) ) 

CALL  VADD(NR,TEMP, A< JJ) ,A(KK) ) 
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6b  KK=KK+NDA 

FAC-DOT <  J ,  U  ( JJ ) ,U< JJ>) 

MR=MR-1 
GO  TO  75 
70  S(J)=1.0 

KK=  1 

DO  72  K=1,JM1 

IF  <  S ( K ) .EQ.l  .  )  GO  TO  72 

TEMP=-DOT(NR,A< JJ) ,A(KK> ) 

CALL  OADD ( K , TEMP  , U ( J J ) , U ( KK ) > 

72  KK=KK+NDA 
75  FAC=1 ./SORT (FAC) 

DO  80  1= JJ , JRM 
SO  A ( I ) =A ( I ) *FAC 

DO  85  I = J J  , JCM 
85  U(I)=U(I)*FAC 
1.00  J J= J J+NDA 

IF  < MR. EQ.NR.OR.MR.EQ.NC)  GO  TO  120 
IF  (MT.NE.O)  WRITE(KOUT, 110)NR,NC,MR 
110  FORMAT ( 13 , 1HX , 1 2 , 8H  M  RANK  ,  12  > 

120  NEND=NC*NDA 
JJ=1 

DO  135  J*1,NC 
DO  125  1=1, NR 
II»I-J 
S  < I ) =0 . 

DO  125  KK- JJ  , NEND , NDA 
125  3< I)=S< I ) +A( II+KK)#U(KK) 

II  =  J 

DO  130  1=1, NR 
U(II)=S(I) 

130  I I =1 I +NDA 
135  J J=J J+NDA 1 
RETURN 
END 

C 

C 

C  HEADER  (IS  A  SUBROUTINE  SINCE  I  GOT  TIRED  OF  LOOKING  AT  IT) 

C 

^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  i L»  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  yU  ^  ^  ^  ■  L.  ^  ^  ^  O, 

^  *  *  *  *  *  *  *  *  *  ®  *  *  *  *  ^  *  *  *  *  ^  ^  ^  ^  ^  ®  ^  ^  ^  ^  ®  ^  ^  ^  ^  *  *  *  *  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ®  ^  ^  ^  ^  ^  ^ 

c 

SUBROUTINE  HEADER 

C 

PRINT' <////>' 

p,  p)  J  ^  *p  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  ^  y^  ^  y^  y^  y^ 

PRINT*,'  *****  *****  ' 

PRINT*,'  *****  DIRECT  OUTPUT  FEEDBACK  *****  ' 
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no  oonnnnn 


PRINT*,'  *****  ***** 
PRINT*,'  *****  ONE  CONTROLLER  VERSION  ***** 
PRINT*,'  *****  ***** 
PRINT*,'  *****  NASA  GROUND  TEST  VEHICLE  ***** 
PRINT*,'  *****  PROGRAM  ACONE  ***** 


PRINT*, '  *********************************************** 

PRINT' <//> ' 

PRINT*,'  THIS  PROGRAM  GENERATES  A  SOLUTION', 

*  '  SUPPRESSING  CONTROL  AND  OBSERVATION  SPILLOVER' 

PRINT*,'  FOR  ONE  SET  OF  SUPPRESSED  MODES,' 

PRINT*,'  RESIDUALS  ARE  ALLOWED  ' 

PRINT*,'  ' 

RETURN 

END 

)>)))))))))))))))))))))))))))))>))))))))))) 

SUBROUTINE  INTEG 

))))))))))>)>>>>>>>))>)>)>))>>))>>>>>>)>>>> 

SUBROUTINE  I NTEG ( N , A , C , S , T ) 

S= INTEGRAL  EA*C*EA  FROM  0  TO  T 
C  IS  DESTROYED 
DIMENSION  A< 1 ) ,C( 1 ) ,S< 1) 

COMMON/MAIN1/  NDIM , NDIM1 ,  X(l) 

COMMON/MAI NB/NCOL,NCOLl 
C0MM0N/MAIN2/CGEF  < 100) 

NN=N*NCOL 
NM1=N- 1 
IND=0 

ANORM~XNORM(N,A) 

DT~T 

IF  ( ANORM*ABS(DT) .LE.0,5)  GO  TO  10 
DT --DT/2 , 

IND=IND+1 
GO  TO  5 

10  DO  15  I - 1 , NN , NCOL 
J-IPNMl 
DO  15  JJ=I,J 
15  S<  JJ)--DT*C<  JJ) 

T 1 =DT**2/2 , 

DO  25  IT~3, 15 

CALL  MMUL (A,C,N,N,N,X) 

DO  20  1  =  1,  N 
II-( 1-1 )*NCOL 
DO  20  JJ=I , NN , NCOL 
11=11+1 


C  <  J J ) - ( X  <  J J ) +X  < 1 1 ) ) #T 1 
20  3  (  JJ  >  =S ( J J ) +C ( J J ) 

25  T 1=DT /FLOAT ( IT ) 

IF  (IND.EQ.O)  GO  TO  100 
COEF ( 1 1 ) =1 .0 
DO  30  1=1,10 
11  =  11-1 

30  COEF (II) =DT#CQEF ( 1 1  + 1 ) /FLOAT  < I ) 

11  =  1 

DO  40  1=1 ,NN,NCOL 
J=I+NM1 
DO  35  JJ=I,J 
35  X< JJ)=A< JJ)*C0EF<1) 

X(II)=X(II) +COEF ( 2 ) 

40  I 1=1 I+NC0L1 
DO  55  L=3 ,11 
CALL  MMUL<A,X,N,N,N,C> 

11  =  1 

T 1 =COEF ( L ) 

DO  55  1=1, NN , NCOL 
J--I+NM1 
DO  50  ,JJ=I,J 
50  X(JJ)=C(JJ) 

X ( I I ) =X ( I I ) +T1 
55  11=11 +NCOL 1 

C  X=EXP<A*DT) 

L=0 

60  L=L+1 

CALL  MMUL (X,S,N,N,N,C) 

11  =  1 

DO  90  1  =  1, N 
J=II 

IF  (I.EQ.l)  GO  TO  75 
DO  70  JJ=I , I I ,NCOL 
S( JJ)=S< J) 

70  J=J+1 

75  DO  85  JJ=I , N 
KK  =  J  J 

DO  80  K=I ,NN , NCOL 
S( J)=S( J)+C(K)»X<KK) 

80  KK=KK+NCOL 
85  J= J+NCOL 

DO  87  JJ=I ,NN,NCOL 
87  C(JJ)=X(JJ) 

90  I I = I I+NCOL 

IF  (L.EQ.IND)  GO  TO  100 
CALL  MMUL <C,C,N,N,N,X) 

GO  TO  60 


75 


n  o  oo  o  n  n 


100  CONTINUE 
RETURN 
END 


SUBROUTINE  L0S<E, 

C2=C2/FL0AT<L) 

C  <  L ) =C2 
D<L+1 )=0. 

I I =N+1 -L 
E  <  1 1 )  =U 
11  =  1 

DO  35  I=1,NN,N 
IL.=  I+Nhl 
DO  30  J=I,IL 
30  X ( J) =EA ( J ) 

X(II)=X(II) +C2 
35  1 1  =  1 1+NF'l 

IF  (L.EQ.N)  GO  TO  40 
CALL  MMUL1(X,A,N,N,N,EA> 

W=W#T1 /FLOAT  <  L ) 

L  =  L  +  1 
GO  TO  20 
40  CONTINUE 

C******  CAN  CHECK  X  0  FOR  ACCURACY 
J=N+63 
DO  50  L=N  t J 
DO  45  K=1 r N 

D<K)=(D<K+1 )-W*C<K) >*T1/FL0AT(L> 
45  E<K)=E<K)+D<K) 

50  U=D ( 1 ) 

11  =  1 

DO  60  I=1,NN,N 
IL=I+NM1 
DO  55  J=I,IL 
55  EA(J)=E<1)#A<J> 
EA(II)=EA(II)+E<2) 

60  I 1=1 I+NP1 

IF  (N.EQ.2)  GO  TO  85 
DO  80  L  =  3 1 N 

CALL  MMUL1<EA,A,N,N,N,X> 

11  =  1 
C2=E ( L ) 


non  nnnnnoon  onnnn 


DO  75  1  =  1 ,  NN  , N 
IL-I+NM1 
DO  70  J=I,IL 
70  EA  <  J ) =X ( J ) 

EA(II)=EA<IIHC2 
75  II=II+NP1 
80  CONTINUE 

THIS  NEXT  SECTION  MULTIPLIES 

THE  MATRIX  EXPONENT  BY  ITSELF  UNTIL  IT  !  ■:  AT 

THE  RIGHT  TIME  STEP 

85  IF  < IND .EO.O)  RETURN 
DO  100  L=1 , IND 
DO  90  I=1,NN,N 
IL=I+NM1 
DO  90  J=I,IL 
90  X ( J ) =EA ( J ) 

100  CALL  MMUL(X,X,N,N,N,EA) 

RETURN 

END 


)  )  >  > 


MLINEQ 


SUBROUT I NE  ML I NEQ <  N , A , C , X , TOL , I ER ) 
SOLOES  A'X+XA+C=0 
A  AND  X  CAN  BE  IN  SAME  LOCATION 
ANSWER  RETURNED  IN  C  AND  X 
DIMENSION  A< 1 ) ,C< 1 ) ,X< 1 ) 
CQMMON/MAINB/  NCOL ,  NC0L1 
C0MM0N/MAIN3/F  < 1 ) 

AD0=T0L#1 ,E-6 
DT=  *  5 
DT1=0 » 

NN=N#NC0L 

DO  5  11=1 ,NN,NC0L1 
5  DT1=DT1-A< II ) 

DTl=riTl/N 

IF  (DTI. GT. 4.0)  DT=DT*4.0/DT1 
11  =  1 

DO  20  1=1, N 
DO  15  JJ=I , NN,NCOL 
15  X( JJ)=DT*A( JJ) 


it  r* vivtvt  ■.-» 


sZ 


■ 


M 

& 
r.  * 
r*; 

f  * 

t  .“• 

i 


ip. 

.*  * 

i 


►  - 


X(II)=X<II)-,5 
20  II=II+NCQL1 

CALL  GMINV(N,N,X,F,MR,0) 

IER  =  4 

IF  (MR.NE.N)  RETURN 
CALL  MMUL(CfF,N,N,N,X) 

C  INITIALIZATION  OF  X,F 
1  =  1 

DO  40  II=1,NN,NC0L 
J=II 

IF  < I . EQ » 1 )  GO  TO  30 
DO  25  JJ=I,II,NCQL 
C< J)=C( JJ) 

25  J=J  +  1 
30  ID=J 

DO  35  J J=I I r  NN »  NCOL 
C  <  J ) =DT*DOT ( N , F ( 1 1 ) , X ( J J ) ) 

35  J=J  +  1 

F(ID)=F<ID>+1,0 
40  1=1+1 

DO  90  IT=1,20 
NEZ=0 

CALL  MMUL(C,F,N,N,N,X) 

1  =  1 
11  =  1 
J=1 

GO  TO  70 
AO  J=II 

DO  65  JJ=I,II,NCOL 
C< J)=C< JJ) 

65  J=J+1 

70  I D= J 

DT1=C( J) 

DO  75  JJ=II ,NN,NCOL 
C ( J ) =C  <  J ) +DOT  <N,F<II),X<JJ>> 

75  J=J+1 

J-J-l 

DO  80  JJ=II. J 
90  X ( JJ ) =F ( JJ) 

IF  (ABS(C(ID)) .GT.1.E150)  GO  TO  95 
IF  (ABS<C< ID) -DTI ) . LT . ( ADV+TOL#ABS ( C< IP U ) ) 
1  =  1  +  1 

11  =  1 I+NCOL 

IF  ( I » LE  * N )  GO  TO  60 
IF  (NEZ.EG.N)  GO  TO  150 
CALL  MMUL (XfXrNfNfNfF) 

90  CONTINUE 
95  IER=1 


NEZ=NEZ+1 


y. 
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RETURN 

150  CONTINUE 
NM1=N-1 

DO  155  I=1,NN,NC0L 
I I^I+NMl 
DO  155  JJ=I,II 
155  X<JJ)^C(JJ) 

IER=0 

RETURN 

END 

)>)>>>)>))>)>>))>>>)>)>))>)>>>>>>)>>>>>>>))>>>)}>)))> 

MHUL 


)>))))))))))))))))))))))))))))))))))))))))> 

SUBROUT I NE  MMUL <  X , Y , N 1 , N2 , N3 , Z  > 

COMMON  /MAINB/  NCOL 

D I MENS I ON  X(NC0L,1) , Y < NCOL , 1 > , Z < NCOL , 1 > 

DO  3  ,.j  - 1  f  N3 
DO  2  I - 1 F  N 1 
S *0 . 

DO  1  K-l f N2 

1  S=3TX< I ,N)*Y(N, J> 

2  Z(I,J)=S 

3  CONTINUE 
END 


-MMUL.1- 


0909090909090909) ( ) ( ) 90 >  < 0 > 9 ) 909 > 9 > 9 > 9 > (09090909 

SIJBRGO  TO  40 
CALL  GMINU(N,N,Z,F,MRf0> 

CALL  MMUL<S,X,N,N,N,Z> 

DO  30  I--=l  ,NN,NCOL 
1 1  =  I  +  NM1 
DO  30  J=I,II 
30  Z( J)=A( J)-Z( J) 

IF  (MR.NE.N)  WRITE ( NOUT y 35 ) MR 
35  FORMAT  ( 27H0RICCATI  SOLN  IS  PSD— RANK  ,13) 

GO  TO  65 
40  CONTINUE 

WRITE ( KOIJT » 45 )  MAXIT 


onnonno  no 


45  FORMAT (27H0RICCAT I  NON-CONVERGENT  IN  ,I2,1.1H 
GO  TO  60 

50  UR  I TE  ( KOUT ,  55 )  IT ,  T 1. 


55  FORMAT (30H0RICCATI  BLOW-UP  AT  ITERATION  , I2,12H 
60  IER=1 
65  RETURN 

200  IF  ( IND . EQ . 2 )  GO  TO  250 
IF  (COUNT. GE. 10.  )  RETURN 
T1=T1/(2.**C0UNT) 

IND=2 
GO  TO  300 

250  T1=T1*(2.**C0UNT) 

INEi-=  1 
GO  TO  300 
END 

>>)))>>>)>>)>>)>>)})>)>)>)>))))>)>)>))>>))>>>>>> 

PRNT 

))>>)))>>)})))))}))})>))))))))>>>>}>>)>  ) )  )  ) 


SUBROUTINE  PRNT < MAT , N ,M > 

COMMON/MA 1NB/NC0L 

COMMON/PRO/Q 

RF MAT  (  NCOL  ,  NCOL  ) 

INTEGER  N,I,J,K,M,Q 
IF  (Q.GT.l)  THEN 
RETURN 
END  IF 

PRINT*,'  ' 

IF  (M.GT.12)  GOTO  2 
DO  1  1=1 , N 

PRINT' (  IX, 12F 10.4) ' , <  MAT ( I , J ) , J=1 ,M) 

CONTINUE 

GOTO  10 

CONTINUE 

IF  (M.GT.NCOL)  THEN 
CALL.  PRNTXL  ( MAT  ,  N ,  M ) 

RETURN 
END  IF 

DO  3  1=1, N 

PRINT' ( IX, 12F10.4) ' , ( MAT ( I , J ) , J=1 , 12) 

CONTINUE 

PRINT' (//) ' 

IF  (M.LE.24)  THEN 


n  n  n 


DO  4  1=1, N 

PRINT' < 1X,12F10,4) ' , < MAT  < I , J ) , J=13,M> 

4  CONTINUE 
PRINT' (///)  ' 

RETURN 

END  IF 

IF  ( M ♦ LE ♦ 36)  THEN 
DO  5  I=1,N 

PRINT' (IX, 12F10.4) ' , <  HAT < I , J ) , J=13,24) 

5  CONTINUE 
PRINT' (//)  ' 

DO  6  1=1, N 

PRINT' (IX, 12F1 0.4) ' , ( MAT ( I , J) , J=24 ,M) 

6  CONTINUE 
ENDIF 

IF  (M.GT.36)  THEN 

PRINT*,'  MATRIX  PRINTING  HAS  BEEN  TRUNCATED 
ENDIF 

10  PRINT' <///>' 

RETURN 

END 


))))))))))) 


>))))))))))))))) 


C  PRNTK 

C  PRINTS  SMALL  MATRICES  WIDER 
C 

C  I))))))))))))))))))))))))))))))))))))))))))))) 
C 

SUBROUT  I NE  PRNTK ( MAT , N , M ) 

COMMON/MAINB/NCOL 
COMMON/PRO/Q 
REAL  MAT ( NCOL , NCOL ) 

INTEGER  N,I,J,K,M,Q 
IF  (Q.GT.l)  THEN 
RETURN 
ENDIF 

PRINT*,'  ' 

IF  (M.GT.10)  THEN 
CALL  PRNT ( MAT , N , M ) 

RETURN 

ENDIF 

DO  1  1  =  1, N 

PRINT'  (  IX,  .1  OF  15.4)  '  ,  <  MAT  ( I  ,  J  )  ,  J=1  ,M) 

1  CONTINUE 

PRINT' (///)  ' 

RETURN 


C  )>>))>>>>>)}>)>>>>>>>>>)>>>>>>)>)>)>)>}))))•<>))> 

c 

C  PRNTXL 

C 

C  )))))))))))))))>))))>)))))>)))))))))))))))))))) 

C 

SUBROUTINE  PRNTXL ( MAT , N , M ) 

COMMON/MAINB/NCOL 
COMMON/PRO/Q 
COMMON/ MA I NA/NDA 
REAL  MAT (NBA, NBA) 

INTEGER  I,J,N,L,M,N,Q 
IF  (Q.GT.l)  THEN 
RETURN 
ENBIF 

PRINT*,'  ' 

BO  1  L-l , M , 12 
K  =  I.  +  11 

IF  (M-L.LT.ll)  K  -  M 
BO  2  1=1, N 

PRINT' (IX, 12F10.5) ' , ( MAT ( I , J ) , J=L , K ) 

2  CONTINUE 

PRINT' (//) ' 

I  CONTINUE 

PRINT' <///)' 

RETURN 

ENB 

C 

t:  ))<><)<><)<><)<)<)<)<)<)<)<)<)<)<)<)<)<><)<><>< 

c 

r;  rfmph.il  this  subroutine  reforms 

C  THE  LINE  OF  SIGHT  PHI  MATRIX 

:  <><><)<><><><><)<)<>()<)()(><><) o ()<)<)()<)(  < 

c 

SUBROUTINE  RFMPHIL < PHIL , RPHIL , NCI , NC2 , NC3 , NR , IC1 , IC2, IC3, IF 
*PETA , CPHIL , MM ) 

COMMON/MAIN A/NBA 
COMMON/MAINB/NCOL 
REAL  PHIL.  ( 3  ,  NCOL 


DIMENSION  X ( 1 ) , A < 1 ) 
COMMON/MAINB/NCOL 
JS=(K-1 )*NCOL*M 
JEND=M*NCOL 
GO  TO  ( .10,30,50,70)  ,  I 
10  BO  20  II=1,N 

BO  20  J.J=I  I ,  JENB ,  NCOL 
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20 


X< JJ)=A< JJ+JS) 

RETURN 
30  DO  40  11  =  1  » N 

KK=< I I - 1 >*NCOL 
DO  40  JJ=1 , M 
LL= ( JJ-1 ) ♦NCOL  +  I I 
40  X  ( NK  + JJ )  =A  ( LL+JS ) 

RETURN 
50  KK=0 

DO  60  11  =  1, JEND  , NCOL 
L.L  =  I  I+N-l 
DO  60  JJ=I  I ,  l.L 
KK=KK+ 1 

60  X(KK)=A< JJ+JS) 

RETURN 

70  KK=M*N+1 

DO  80  11  =  1  ,  M 
I.L=(M-II  )*NC0L+1 
DO  80  I J  = 1 , N 
KN=KK-1 
J J=LL+N-I J 
30  0 ( JJ+JS ) =X ( KK  ) 

RETURN 

END 

C 

C 

C  TIMEX 
C 

C  $#**$***$#$***$#*#$#$***$#$#$#$#$#$*$#$*$*$#$#$**#$# 

C 

SUBROUTINE  TIMEX  <  STM , MM , DT , XO , PUT , TMAX , X 1 , XC , PHILOS , 
*IC1 ,  IC2,IC3,NC.l,  NSUP ,  NRES  > 

C0MMQN/MAINA/NDA,NDA1 
COMMON/MAI NB/NCOL, NCOL 1 

REAL  DT , TMAX , ABT , STM( NDA , NDA ) , XO ( NDA ) , XI ( NDA ) ,XC(NDA) 

C  OK, SO  THIS  NEXT  LINE  IS  NOT  THE  ELEGANT  WAY  OF  DOING 

C  THINGS  IS  YOUR  THESIS  DUE  TOMORROW? 

C  E  IS  MY  OUTPUT  VECTOR  WHEN  YOU  CALL  LOS<  )  THE  NUMBER  IS 
C  THE  SIZE  OF  YOUR  OUTPUT  VECTOR.  MINE  IS  THREE 
C 

REAL  E ( 3 ) 

INTEGER  I , J,MM , PDT ,EN 
EN  =  3 
ABi  0. 

C  MAKE  A  COPY  OF  THE  I.C.  VECTOR  AND  USE  THE  COPY. 

DO  10  1=1, MM 
XC<I>  =  XO(I) 


nnonnnoon 


10  CONTINUE 
PRINT' (//>  ' 

PRINT*,'  I  NIT  STATE  <  TIME  =  0.0)' 

PRINT' (IX ,  10E12  ♦  4 ) '  ,  (XOd  ) » 1=1 , MM ) 

CALL  LOS (E, PHILOS, EK, XO, IC1 , I C2,IC3, NCI, NSUP,NRES> 
PRINT*, AST ,E<1) ,E<2> ,E<3) 

C  CHECK  THAT  MAX  TIME,  ' TMAX ' ,  IS  NOT  REACHED  AND  CONTINUE  TO 
C  PROPAGATE  THE  STATE  VECTOR  AND  PRT/PLT  IT  EVERY  ' PDT '  TIME 
C  INCREMENTS 

20  DO  40  1=1, PDT 

CALL  VMULFF ( STM , XC , MM ,MM, 1 , NDA , NDA , XI , NDA , IER > 

DO  50  J=1 , MM 
XC<J>  =  X 1 ( J ) 

50  CONTINUE 
40  CONTINUE 

C  INCREMENT  THE  ABSOLUTE  TIME  AND  OUTPUT  THE  STATE  VECTOR 
ABT=ABT+DT*PDT 

CALL  LOS ( E , PHILOS , EK , XC , IC1 , IC2 , IC3 , NC 1 , NSUP , NRES ) 
PRINT* , ABT , E  < 1 ) ,E<2> ,E<3> 

IF(ABT.LE.TMAX)  THEN 
GOTO  20 
END. IF 
RETURN 
END 

*<*<*<*<*<*<*<*<*<*<*<*<*<*<*<*<*<*<*<*<#<*< 

-TOUT- 

THIS  SUBROUTINE  PRINTS  THE  TIME  RESPONSE  AS  REPRESENTED 
BY  SENSOR  OUTPUTS 

J*)*)*)*)#)*)*)*)*)*)*)*)*)*)*)*)*)#)*)*)*)*)* 

SUBROUTINE  TOUT (Cl , NCI ,C2, NSUP, C3, NRES, XC,Y> 

C  NEED  COMMON  BLOCKS  FOR  NSEN 
COMMON/M A I NA/NDA,NDA1 
COMMON/M A I NB/NCOL ,NCQL1 
COMMON/ 10 /NSEN , NACT 

REAL  Cl  (NCOL,NCOL)  ,  C2  <  NCOL  ,  NCOL  >  ,C3(NC0L  .  ’HI.) 

REAL  XC(NDA) , Y(NSEN) 

DO  1  1=1, NSEN 
Y ( I ) =0 . 0 
1  CONTINUE 

DO  2  1=1, NSEN 
C 

DO  3  J=1 » NCI 

Y(I)  =  Y(I)  +  Cl < I , NCI + J )  *  XC ( NC 1 + J ) 

3  CONTINUE 


-1  i»'in  »v« ^■^rrr^T  ■  ■  ■  ■  I  ■!■»-■■  t  ' L  * i 1 


B 


I 


► 


B  <• 


C 

C 

c 

c 

c 

c 


c 

c 

c 

c 

C 

c 

c 


K  -  <NC1*2)  +  NSUF' 

DO  A  J  =  1 f NSUP 

Y(I)  =  Y < I )  +  C2(Ir NSUP+ J )  *  XC(K+J) 

A  CONTINUE 

K  =  <NC1*2)  +  <NSUP*2)  +  NRES 
DO  5  J  =  1 f NRES 

Y ( I )  =  Y(I)  +  C3 ( I r NRES+J )  *  XC(K+J) 

5  CONTINUE 
2  CONTINUE 
RETURN 
END 

))))))))))))))))))))))))))))))))))))))))))))))  >)  ) ) 

VADD 

ADDS  VECTORS 

SUBROUTINE  VADD < N , Cl , A f B ) 

DIMENSION  Ad)  ,Bd) 

DO  1  I - 1 f  N 
Ad  )=Ad  )  FCl*Bd  ) 

RETURN 

END 

XNORM1 

)))))))))))>)))))))))))))))))>))>)))))>)))>))>))>)>)* 
FUNCTION  XN0RM1<N,A) 

COMPUTES  AN  APPROXIMATION  TO  NORM  OF  A —  NOT  A  B  ID 
DIMENSION  A<1) 

COMMON/MA INA/  NDA,NDA1 

NN=N*NDA 

C1»0» 

TR  =  A< 1 ) 

IF  (N.EQ. 1 )  GO  TO  20 

r~2 

DO  10  I  I==NDA  1  r  NN  f  NDA 
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ESOLUTION  TEST 


*  %  m  -  t 

V  .*•’.*** 

V.'.V, 

V.v.v 

s.*  •;  %' 

v‘>  v 

i-VV 

yvv 


THIS  FILE  CONTAINS  THE  INFORMATION  NEEDED  TO  RUN  ANY  ACOSS  DERIVITI 
PROGRAM  FOR  THE  NASA  GROUND  TEST  FACILITY  APPARATUS 
THIS  INFO  IS  FOR  THE  ANTENNA  CONFIGURATION  WITH  THE  VCOSS 
MODS.  THIS  INFO  HAS  NOT  EXPERIMENTALLY  VERIFIED,  SO 
IT  REALLY  IS  BORDERLINE  UTILITY 
DOF  1. 

. 1016E-02  -.6448E-01  .34S6E-06  -.7749E-01  -.1A13E-05  -.3119E-03 
.32S1E-04  . 1832E-02  .7944E-.7  -.6870E-03 

-.64S4E-12  1673E-11  .1342E-01  .3390E-03  .3294E-02  .U18E-08 

- . 3998E-04  -.9612E-07  -.2A34E-17  -.2228E-16  -.1828E-0S 

-.2506E-03  -.3083E-01  .4494E-04  .4301E-01  -.3I83E-05  -.I282E-04 


-.6434E-12 
- .3998E-04 
-.2506E-03 
1 .580E-01  . 

DOF  6 
.1A99E-14 
-.97.3E-13 
•7131E-10  . 

-.1144E-24 
.3717E-09  - 

OOF  8 

-.2247E-04 
-.2098E-03 
.7300E-01 
-.3287E-16 
•4904E-01 
DOF  12 
-.1808E-20 
-.2S2SE-10 
-.6743E-14 
.177SE-23 
-.1815E-12 
DOF  13 
.2022E-12 
-.433AE-01 
1036E-04 
.2389E-15 
- . 9389E-03 
DOF  19 
. 2022E-1 2 
- . 7914E-02 
.2073E-01 
1045E-04 
•3079E-08 
-.131AE-01 
- . 9389E-03 
.9A17E-04 
DOF  20 
-.2247E-6 
.7999E-04 
. 3080E-06 
.7300E-01 
-.2034E-03 
-.3932E-03 
.4906E-01 
. 1743E-07 
DOF  21 


•1S31E-08  -.3329E-03 


. 1566E-18 
.6421E-U 
.1931E-13 
•4360E-23 
-.3884E-13 

.2981E-10 
• 1034E-01 
.1733E-04 
.6379E-17 
-.4678E-0S 


-.6380E-21 
. 1848E-15 
.8703E-10  . 

1613E-13 
.4273E-13 

-.2384E-11 
.3080E-06 
• 1633E-01  . 

- .3952E-03 
. 4796E-05 


.1144E-12  .5347E-15  -.4409E-11 

-.1808E-11  -.2993E-20  -.3036E-20 

3480E-1A  -.1245E-11  -.S380E-13 

-.2319E-13  -.3606E-09  .4938E-14 

-.8683E-09  • 3861E- 17  .3233E-11 


.1424E-03 
-.2744E-02 
.3688E-08  - 

-.7406E-03 
-.S493E-01 


.7999E-06  -.A712E-02 

-.4033E-11  -.7391E-11 

•  2034E-03  -.3433E-06 

1122E-00  .9482E-06 

. 1945E-09  .2723E-03 


- . 4804E-20  -.3846E-14  .3321E-16 

-.2343E-14  .1308E-10  .2312E-13 

.3007E-10  .3314E-14  -.7A63E-14 

.2237E-24  -.3212C-10  -.5S33E-10 

-.8602E-09  -.2363E-08  -.2009E-12 


-.9801E-11  -.1212E-14 
.9799E-21  -.1033E-19 

-.7220E-16  -.2210E-09 

-.1079E-I3  -.3227E-12 

.2229E-11  .2203E-13 


.1231E-10  .1634E-08  -.2393E-06  -.7914E-02  -.1077E-05 

-.4058E-03  . 2073E-01  .3743E-04  .1449E-11  -.1436E-10 

.4832E-01  -.3022E-06  . 7033E-06  .5079E-0B  -.2940E-01 

•3531E-16  1526E-01  .6049E-02  -.582BE-06  -.1801E-04 

-.3341E-01  1238E-00  -.9432E-03  .9617E-04  ,9649E-0e- 


.1231E-10 

1077E-03 

.3743E-04 

-.4832E-01 

-.2940E-01 

•6049E-02 

-.3341E-01 

.9449E-06 

.2981E-10 
-.A712E-02 
-.2744E-02 
.1733E-04 
-.34S3E-04 
- • 7  40AE-03 
- . 4A7RE-0S 
.2723E-03 


.  1634E-03  2393E-06 

-.4336E-01  -.4058E-03 

01449E-11  - • 1 636E- 10 

-.3022E-06  .7033E-6 

.2389E-15  .3331E-IA 

-.3328E-06  1 30  IE -04 

-.1238E-00  -.9432E-03 


- . 2334E- 1 l 
-.207SE-03 
— .4030E-11 
.1633E-^1 
-.3287E-14 
-.1122E-00 
,479oE-05 


. 1 424E-03 
.1034E-01 
-.7391E-11 
3488E-0S 
.6179E-17 
.9482E— >6 
-.3  193E-01 


*-  •  .  %  • 

.v: 


\v.-v 

■Vs*: 


.1016E-02  -.6448E-01  .3636E-06  -.7749E-01 

-.1613E-03  -.6712E-02  .3231E-06  . 1832E-02 

.7944E-7  -.6870E-03  -.6434E-12  -.1673E-11 

.1342E-01  .3390E-3  .3294E-02  .1118E-08 

- . 3998E-04  - • 96 12E-07  -.2636E-17  -.2228E-16 

-.1828E-03  - . 2306E-3  -.3085E-01  .4494E-06 

•4301E-01  -.5183E-03  -.1282E-04  .1380E-00 

•  1331E-08  - « 3329E-03 

DOF  46 

• 36E-1  .2323E-6  .2068E-7  -.1611E-2  -.145E-3  .U92E-1 
•3964E-6  -.1630E-2  -.4489E-7  .4108E-3  .6283E-12  .1034E-11 
-1.003E-2  -.2314E-5  -.3803E-2  -.1492C-B  .3338E-4  .1632E-7 
.6333E-17  -.1601E-13  .3843E-6  .3163E-6  .8403E-2  -.7136E-7 
-.3871E-2  .3739E-6  -.2788E-6  .3337E-2  -.2162E-10  -.1B11E-4 
DOF  47 

•1917E-7  - . 189E-6  -.3339E-1  .1741E-6  -.9893E-2  -.1167E-3 
-.3841E-2  -.2243E-6  .1184E-2  . 1936E-3  .7633E-13  -.8048E-12 
-.3331E-6  .1342E-2  .3394E-7  -.1224E-6  -.1133E-8  -.2406E-2 
-.6832E-16  ■ 3606E-17  -.4618E-3  -.4436E-3  -.6023E-7  -.1737E-S 
-.7039E-6  -.31E-2  -.6223E-2  -.4599E-6  .4648E-5  .4393E-7 
DOF  31 

-.18T4E-10  • 2788E-10  .3413E-06  .1086E-5 

-.1195E-01  -.109’E-03  .3139E-01  .2932E-03 

-.9038E-02  -.1604E-04  -.1000E-12  .3021E-11 

-.3376E-03  . 4583E-01  -.3997E-03  .2267E-04 

•8119E-07  -.l’09E-02  - . 1 142E-1S  .3963E-17 

-.1608E-01  -.1193E-01  -.3346E-03  -1022E-03 

-.2721E-04  -.1929E-02  -.7346E-02  .8200E-03 

-.3337E-04  -.236OE-03 


mi 

m 

wi 

$$ 

-\v-- 

■w  "L 


18T4E-10 
-.1193E-01 
-.9038E-02 
-.3376E-03 
•8119E-07 
-.1608E-01 
-.2721E-04 
-.3337E-04 
DOF  32 
-.4123E-06 
.1241E-03 
-.6721E-06 
-.6333E-01 
1063E-02 
-.2876E-03 
2684E-01 
-. 1369E-07 
DOF  53 
.1012E-02 
. 123AE-04 
-.2193E-05 
. 1 773E-02 
-.3978E-04 
• 3301E-06 
-. 30S2E-01 
-.1298E-06 


-.9032E-10 
- . 1203E-01 
.4940E-02 
-.1177E-04 
•8967E-07 
-.3337E-03 
.3270E-03 
-.4686E-03 

-.6448E-01 
*  4233E-02 
-.2403E-01 
-.4366E-04 
•4324E-06 
- . 1948E-03 
.2043E-03 
-.8707E-03 


-.4681E-11 

.4866E-3 

.3317E-11 

.7074E-01 

-.6228E-16 

-.7839E-01 

.9339E-03 


. 1949E-04 
-.1343E-01 
. 1600E-10 
.2963E-07 
-.3497E-16 
.6143E-06 
-.3133E-01 


*  3A33E-OA 
4414E-04 
-.9792E-11 
. 1242E-02 
•4623E-18 
-.1186E-01 
-.1330E-04 


.3343E-01 
-.3933E-02 
.4148E-19 
- i B209E -08 
.4010E-13 
.  1826E-07 
.6633E-02 


THESE  ARE  THE  COORDINATES  THAT  DESCRIBE  THE 

POSITION  OF  THE  ANTENNA,  66  67  63  ARE  X,Y,Z  RESPECTIVELY. 

69,70,71  ARE  THETA  Y,Y,Z  RESPECTIVELY. 

DOF  66 

. 3803E- 1  -.1309  .7S22E-6  .107SE-0  .3916E-4  -.1131 
-.3909E-4  - . 227E- 1  -.1424E-5  .4304E-1  -.16»2E-10  . 1043E-9 
- . 2386E -1  -.11 13E-3  .2308E-1  -.2239E-'  -.513CE-3  .3131E-5 
- . 5697E - 1 7  . 3695E-13  .2777F-4  .l'!9E-4  -.3221E-1  . 33 1 7E  -4 
-.9199E-1  .  1 733E-4  .2704E-4  .U44E-1  -.2594E-'  -,2cj48E-2 


89 


OOP  47 

.1924E-7  - . 1893E-4  -.334E-1  -.2819E-S  .1204  .1397E-4 
-.1444E-1  - . 3427E-5  -.1934E-2  .1242E-4  -.1  48E-U  -  •  '-11 

-.3148E-3  . 932E-2  .1237E-5  -.S11E-4  -.304E-7  -.1042E-1 
• 1231E-14  .4318E-17  .8381E-2  -.1044E-1  . 2998E-4  .7942E-3 
•4333E-S  .1807E-2  .4231E-2  -.1401E-3  -.3478E-3  .4871E-6 
OOP  48 

.3782E-10  -.3734E-10  -.4914E-4  -.220BE-3  .2437E-1  .2233E-3 
-.1101E-0  -.7387E-3  .1322E-1  .3242E-4  -.2334E-12  -.3092E-11 
.1203E-4  - • 1035E-0  .1249E-4  .4S33E-3  - . 2444E-4  .4749E-2 
.1404E-13  -.8392E-17  .1173E-0  .1309  -.4340E-4  .3274E-3 
-.9317E-4  -.3474E--  -.9442E-2  .93S3E-3  -.2742E-3  -.2338E-3 
DOF  49 

•• 1 849E-10  .2849E-10  .3401E-4  .1088E-3  -.1202E-1  -.1102E-3 


•3383E-1  .4209E-3  -.3784E-2  -.1334E-4  -.2334E-12  -.S092E-U 
-.4138E-3  .3447E-1  -.7249E-3  -.27E-4  .1411E-4  -.3344E-2 
* • 4384E-14  .3339E-17  -.7373E-1  -.8473E-1  . 4843E-4  -.324E-3 
•7978E-4  . 4847E-2  .1031E-1  -.108E-4  .7134E-4  .2833E-5 
DOF  70 

-.4137E-4  -.7492E-10  -.4397E-11  .904E-3  .1222E-5  -.1199E-1 
.3748E-3  -.2387E-1  -.3133E-4  .4809E-2  .4809E-11  .2244E-10 
-.4343E-1  -.1233E-4  .7001E-1  .212BE-7  -.9B99E-3  .1089E-4 
-.4974E-14  -.3138E-14  -.2498E-5  -.32S3E-3  -.7931E-1  .4228E-4 
-.2793E-1  . S497E-5  .9984E-3  -.3302E-1  .1134E-7  . 184E-3 
DOF  71 

.1012E-2  -.4448E-1  • 3433E-4  .S334E-1  .1274E-4  .3734E-2 
-.4843E-4  -.3127E-2  .1833E-3  .1741E-1  -.4738E-11  .3379E-10 
-.8287E-3  -.3328E-4  .7B43E-2  - . 1774E-7  -.1308E-3  .1434E-5 
- . 53S3E-17  .4103E-13  .1343E-4  .7324E-3  -.4993E-1  .4237E-4 
-.1234  .2133E-4  . 3231E-4  .2417E-1  .S812E-7  -.178E-2 
DOF  28 


.036 

. 232E-4 

.2044E-7 

-.9448E-3 

.3144E-3  - 

>  4437E-1 

-• 1 103E-4 

■3134E-1 

.  M39E-3 

-. 1311E-1 

- . 1942E-10  - . 3401E-10 

*  2892E0 

. 6 1 04E-4 

-.2284E0 

- . 1008E-4 

• 3403E-2 

.2310E-3 

.2174E-15 

-.2280E-13 

-.2037E-4 

-. 1444E-4 

.5333E-1  - 

.4990E-4 

“ » 1624E0 

. 1923E-4 

. 4092E-3 

.4043E-1 

-.3818E-8 

.3433E-3 

OOF  29 

. 19X0E-7 

-.1B91E-4 

3339E- 1 

. 1 173E-5 

. 4681E-3 

.4230E-5 

. 1770E00 

.1144E-4 

-♦3944E-1 

-.9441E-4 

-33422E-11 

. 3831E-10 

.6400E-3 

.1739E-1 

-.1446E-4 

. 3749E-4 

. 2987E-6 

.  2029E00 

- . 1331E-14 

-.2143E-13 

1338E00 

. 2784E00 

. 1037E-4 

.2738E-3 

» 3071E-4 

- . 4109E- 1 

- . 8299E-1 

-. 1220E-4 

. 1902E-3 

.  2101E-3 

DOF  40 

. 3600E-1 

.2318E-6 

•  2063E-7 

- • 7745E-3 

. 1264E-4  - 

.  1137E00 

.2233E-5 

- . 2481E-2 

- • 61 48E-6 

. 2038E-2 

- .3102E-1 l 

.1783F  10 

. 1423E-1 

- . 1943E-6 

-.2963E-1 

1089E-7 

. 3817E-3 

•5887E-4 

. 4O07E-16 

-.4293E-13 

.  1034E-4 

.1173E-4 

.7184E-1 

.  4337E-4 

*  1 161E0 

-.  1437E-4 

- » 1372E-4 

- . 9338E-2 

* • 1 668E- 7 

.  4678E-3 

DOF  41 

• 1923E-7 

- . 1893E-6 

- • 3340E-1 

- . 2136E-3 

•  1 133E00 

.  1330E-4 

.  1301E-1 

1827E-3 

7468E-2 

. 2840E-5 

-.1129E-11 

- . 3088E-1 1 

-.6446E-3 

.3726E-1 

- .  1 174E-3 

- • 3728E-4 

.1890E-7  - 

.  U44E-1 

. t 131E-14 

.6912E-17 

-.1240E-2 

- . i 768E- l 

- . 1 790E-3  - 

.5432E-4 

1218E-4 

.6229E-3 

-.2494E-3 

• 3622E-3 

-  *  5890E-3  - 

.  9380E-4 

OOF  48 

. 3600E-01 

.  2318E-06 

.  2063E-07 

- . 7626E-03 

. 1 339E-04 

-1 . 21  IE-' 
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.5203E-05  - 

.1373E-01 

102E-03  . 

5071E-02  . 

1 408E-12  .2 

:~aie-io 

- . 259 IE -01 

-.7373E-05 

. 1 350E-01 

. 7201E-0B 

- • 266oE-03 

.6434E-06 

.2874E-17  ~ 

. 4506E-15 

. 3582E-05  . 

3487E-05 

2391E-01 

5902E-07 

.9976E-01  ' 

.  1 336E-04 

-.7999E-0S 

-.2857E-01 

- . 2626E-07 

.3B19E-03 

DOF  49 

.1926E-07  - 

» 1B93E-06 

- • 3340E-01 

-.2819E-05 

1 .206E-01 

. 1397E-04 

1646E-01 

-.3627E-05 

1934E-02 

.1262E-04 

- • 1068E-1 1 

- . 6152E- 1 1 

-  4  3 168E-05 

.9320E-02 

.  1257E-05  - 

.5111E-04 

-.3040E-07 

-.  1042E-01 

.  1 251E- 14  . 

43X7E-17 

.  8580E-02  -. 

1044E-01  . 

299QE-06  .7960E-03 

.  4T/34E-03  . 

180 6E-02 

.  4250E-02 

1400E-05  - 

. 5684E-03  . 

407OE-O6 

DOF  30 
.2810E-14  - 

* 1726E- 1 3 

.4308E-14  - 

.  1202E-10 

. 6696E-06  . 

6404E-10 

1837E-04 

-.3464E-08 

-. 1030E-04 

.  1989E-08 

-.5700E-14 

.  2524E-13 

.1548E-08  - 

•8275E-05 

- .  8377E-07 

. 158SE-05 

.2010E-08  - 

• .  1212E-02 

.4604E-17 

t 1 339E" 17 

.641 IE-02 

.5109E-02 

-.2248E-05 

.  1Q86E-04 

- • 3723E-05 

- . 2301E-03 

- • 5525E-03 

. 5904E-06 

-.3228E-04 

-  9 1595E-06 

NATURAL  FREQUENCIES  RAD/SEC 

0  0 

0  .3601389  .8209141 

.8404761 

2.11471 

2.665896 

2.74208 

2.835489 

2.950763 

2.954657 

3.604164 

3.903844 

6.464519 

6.526102 

6.526102 

7.903797 

8.15414  8. 

.16517 

8.761278 

1  9.433749  11.04083 

11.24278 

13.24458 

13.7360 

16.7332 

17.23949 

18.06101 

18.06101 

NATURAL  FREQUENCIES  SQUARED. 

0.0  0.0  0.0  . 1297  . 6739  .706*  *.*72  7.107  7.519  8.040  8.707 

8.730  12.99  15.24  41.79  42.59  42.59  62.47  66.49  66.67  76.76 

89.09  121.9  126.4  175.4188.7  280.0  297.2  326.2  326.2 

•1016E-02  -.&443E-01  . 3656E-06  -.7749E-01  -.1615E-05  -.3119E-03 
-.1615E-05  -.3U9E-03  .3251E-06  . 1852E-02  .7944E-.7  -.6870E-03 

-.a454E-12  -.1675E-11  . 1342E-01  .3390E-05  .3294E-02  .1118E-0B 

-.3998E-04  -.9612E-07  -.2636E-17  -.2228E-16  -.1B28E-05 

-.2506E-05  -.5085E-01  .4494E-06 
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Abstract 


Direct  output  feedback  control  using  one  or  two  controllers  is  applied  to 
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used  to  design  the  gain  matrices,  and  full  state  feedback  is  approximated  by 
use  of  generalized  inverses  of  the  observation  matrices.  Spillover  is 
eliminated  through  the  use  of  transformation  matrices. 

The  structure  is  shown  to  be  controllable  with  this  method.  Alternative 
sensor  placement  is  explored,  and  found  to  cause  improvement  in 
performance.  The  torsion  modes  are  found  to  be  particularly  important  to 
the  performance  of  the  structure,  but  need  more  sensing  and  actuation  to  be 
adequately  controlled.  Two  controller  systems  require  more  sensors  and 
actuators  than  available  to  achieve  acceptable  performance  with  this 
structure. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


